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The resonant structure of Jupiter's trojan asteroids-II. 
What happens for different configurations of the planetary 
system. 
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ABSTRACT 

In a previous paper, we have found that the resonance structure of the present Jupiter 
Trojan swarms could be split up into four different families of resonances. Here, in a 
first step, we generalize these families in order to describe the resonances occurring 
in Trojan swarms embedded in a generic planetary system. The location of these 
families changes under a modification of the fundamental frequencies of the planets 
and we show how the resonant structure would evolve during a planetary migration. 
We present a general method, based on the knowledge of the fundamental frequencies 
of the planets and on those that can be reached by the Trojans, which makes it possible 
to predict and localize the main events arising in the swarms during migration. In 
particular, we show how the size and stability of the Trojan swarms are affected by 
the modification of the frequencies of the planets. Finally, we use this method to study 
the global dynamics of the Jovian Tr ojan swarms when Satu rn migrates outwards. 
Besides the two resonances found by iMorbidelli et alj (|2005f ) which could have led 
to the capture of the current population just after the crossing of the 2:1 orbital 
resonance, we also point out several sequences of chaotic events that can influence the 
Trojan population. 

Key words: celestial mechanics - minor planets, asteroids - Solar system: general. 



1 INTRODUCTION 

The discovery of Achilles by Wolf in 1906, and of four 
other Jovian Trojans the next year, gave a new impulse 
to the study of the triangular configurations of the three- 
body problem, whose existence was shown by Lagrange 
in 1772. An important problem was to establish the ex- 
istence of a stable area surrounding the triangular equi- 
librium points La and L5 associated to the Sun and 
Jupiter. From a mathematical point of view, the appli- 
cation of K.A.M. theory to the planar and circular re- 
stricted three-body problem, gives a result of confinement 
between two invariant tori, ensuring the sta bility in the 
neigh b orhood of L4 and L5 for an infinite t i me (|Leontovitchl 
19621: iDeprit fc Deprit-Bartholomel 1 19671 : iMarkeevI 1 19721 : 
Meyer fc Schmidd ll986l . In the case of the spatial restricted 
three-body problem, where K.A. M. theory does not en sure 
infinite time stability anymore, iBenettin et all (Il998l) ap- 
plied an extension of the Nekhoroshev theorem |l977j) to 
quasi-convex Hamiltonians in order to prove exponentially 
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long-time stability. But these two complementary theories 
do not give any information about the size of the stable 
region surrounding the triangular equilibrium points. In or- 
der to get estimates of this width, several authors developed 
Nekhoroshev-like estimates based on normalization up to 
an op t imal order ijGiorgilli et al ] |l989l; ICelletti fc Giorgiul 
I l99ll : iGiorgilli fc Skokosl [19971 ISkokos fc Dokoumetzidij 
I200H ). More recent results, based on improvements of 
these methods, can be found in ([Gabern et al . 2005; 
lEfthvmiopoulos fc Sandon 120051 ; ILhotka et al l 120081) . An - 
other very interesting result was published by Rabe 1 1967h . 
In this work, based on the change of stability of the family of 
long-period Lyaponov orbits emanating from L4, the author 
gives an approximation of the limit eccentricity of stable 
tadpole orbits with respect to their amplitude of libration. 

Except Gabern and Jorba wh o applied the same meth- 
ods as in i|Giorgilfi et all Il989h to the bicircula r , tri- 
circular and biannula r coh erent problem s llGabernl 120031 ; 
iGabern fc Jorball2004r t, and ILhotka et all (|2008l ) who used 
a more sophisticated method in the elliptic restricted three- 
body problem (RTBP), all of the above mentioned works 
were developed in the framework of the circular RTBP, pla- 
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nar or spatial, which is not a very realistic model for the pur- 
pose of studying the long-term dynamics of the Jovian Tro- 
jans. To overcome this proble m, several authors perform ed 
full numerical integrations. In ll993l . iHolman fe Wisdom! , in- 
tegrating test-particles during 20 Myr, obtained the first 
global stability result concerning the Trojan s of the giant 
planets in our Solar system. The same year, iMilanil devel- 
oped a numerical method providing the proper elements of 
the Trojans of Jupiter. A semianalytical model of deter- 
mi nation of these pr oper e lements has also been proposed 
by Beauge fc Roid in l200ll . Four years later, iLevison et al.l 
i| 19971 ) described the long-time erosion of the Jovian Tro- 
jan swarms without identifying its mechanism (more re- 
sults concerning the long-term beha vior of the Trojans 
are published in iTsiganis et al.l (120051)). A few years later, 

Donesl (|2002T l 



Michtchcnko et al 
Marzari et all (200. 



Nesvornv fc 



and 



20011 ) 

0, showed the existence of two sets of 



unstable structures lying inside the Jovian Trojan swarms 
which could be related to the long-time erosion. The authors 
suggested that one of these sets was connected to the great 
inequality (between Jupiter and Saturn) while the second 
one may have been generated by the commensurability be- 
tween the libration frequency of the co-orbital reson ance and 
the frequency n.j up — 2ns a t- In 20051. iRobutel et ail identified 
these singularities and analyzed their underlying resonances, 
this leading to the decomposition of th e "resonance struc- 
ture" in four families of resonances. In l|Robutel fc Gabernl 
2006), hereafter called Paper I , a new description of this res- 
onant structure was given, and the link between these reso- 
nances and the long-term stability (or instability) of Jovian 
Trojans was established. 

The structures described in Paper I , generated by com- 
mensurabilities between the proper frequencies of the Tro- 
jans and the fundamental frequencies of the planetary sys- 
tem, depend on the planetary configuration. Therefore, a 
small modification of the geometry of the system is able 
to modify the Trojan swarms' resonant structure, chang- 
ing the global dynamics, and, consequently, the stability of 
the co-orbital region. Plan e tary m igration in a planetesi- 
mal disc (see I Gomes et alj {2004) and references therein) 
provides a natural m echanism of evolution of planetary 
systems. According to iMorbidelli et al.l l|2005l ). the Jovian 
co-orbital population is not primordial, but was, instead, 
captured during chaotic events which took place in the 
course of a planetar y migration. The so-called Nice model 
|Tsiganis et al.ll2005l ) predicts that during the migration (in- 
wards for Jupiter, and outwards for Saturn, Uranus and Nep- 
tune), t he Jupiter-Saturn coup le crossed the 2:1 orbital res- 
onance. IMorbidelli et al] l|2005l ) suggested that the Trojans 
were captured in co-orbital motion just after the crossing of 
this resonance, when the dynamics of the Trojan region was 
completely chaotic. The authors show that this global chaos 
arises when the libration frequency of the Trojans is close to 
the combinations 3(nj up — 2ng at ) or 2(nj uv — 2ns a t )- Thi s 
result was recently confirmed by iMarzari fc Scholil l|2007l ) , 
showing that these resonances are the main generators of 
the depletion (and, probably also the capture) of the Trojan 
population. 

The results reported in P aper I are directly relat ed to 
the phenomenon descri bed in IMorbidelli et all (|2005l ) and 
IMarzari fc Scholil l|2007l ). More precisely, the resonances in- 
volved in this process are members of one of the four families 



of resonances presented in Paper I . Consequently, this pre- 
vious paper contains the necessary material to derive tools, 
that make it possible to study the evolution of the Trojans 
resonant structure. This evolution is induced by changing 
the geometry in the considered planetary system. Indeed, 
our goal is not to give an accurate and realistic description 
of the behavior of the Trojan swarms during a migration pro- 
cess, but rather to study what are the planetary configura- 
tions for which the Trojan swarms become globally chaotic. 
Moreover, we want to develop a model that is capable of 
making predictions in a large class of problems, and that is 
not restricted to the study of the Jovian Trojans. Our ap- 
proach is based on three different points. The first one rests 
on the understanding of the Trojans' resonant structure, and 
its decomposition. It will be shown in section [5] that, for 
a generic Trojan swarm, almost all resonances driving its 
global dynamics are members of four different families which 
generalize the families described in Paper I . The resonances 
of all these families involve both fundamental frequencies 
of the Trojans and the basic planetary frequencies. These 
families establish a link between the behavior of the Trojan 
swarm and the geometry of its planetary system. The two 
other ingredients on which the method is based are respec- 
tively the exploration of the planetary frequencies and the 
determination of the frequency domain to which the main 
frequencies of the whole Trojan swarm belong. In section |3j 
our method is applied to the study of the Jovian Trojans per- 
turbed by Saturn, employing the model used in Pape r I but 
also in (|Morbidelli et al]|2005l ; IMarzari fc Schol]||2007l ). This 
application reveals the rich dynamics of the Trojan swarms 
under the action of numerous resonances of different origins, 
which generate a resonance web which cannot be fully un- 
derstood without knowing the fundamental frequencies of 
the considered dynamical system. 



2 THE RESONANCE STRUCTURE OF A 
GENERAL TROJAN SWARM 

2.1 General setting 

Let us consider a general case: the planet harboring the 
Trojan swarms is the p th planet of a given planetary sys- 
tem composed of N planets orbiting a central star. We 
denote by rrij, respectively Oj, the masses and the semi- 
major axes of the planets, and mo the mass of the star. 
The small mass parameter p is defined by the expres- 
sion p = max(mi/M, ••• , mjv/M), where M is the total 
mass of the system. A linear combination of secular fre- 
quencies will be denoted: k g • cr g + l s ■ cr s where k g and 
l s are elements of Z . If we assume that the planetary 
system is stable enough to be considered as quasiperiodic 
for a given time length (see paper I), its fundamental fre- 
quencies are denoted: (ni,--- ,njv), <r g — (<?i,--- ,<?jv) and 
f s = (si, • • • , sjv) El where the rij denote the proper mean 
motions which are of order with respect to the planetary 
masses. The vectors cr g and cr s correspond to the secular 
frequencies of the planetary system (respectively associated 



1 Owing the invariance of the total angular momentum, one of 
the Sj is equal to zero. 
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to the precession of the perihelia and to the precession of 
the ascending nodes) which are both of order /i. 

2.2 Restricted three-body problem 

In a first step, we retain only the gravitational interaction 
of the p'th planet (and of the star) on the Trojans. We 
are thus brought back to the RTBP. This model gives, at 
least for the outer planets of the Solar system, results which 
are quite realistic (this is not the case for inner planets, 
see section I2,3|l . In this model, the triangular equilibrium 
points Z/4 or L5 are well defined, even for an elliptic mo- 
tion of the planet (paper I). Strictly speaking, in an inertial 
frame, the equilateral configurations correspond to periodic 
orbits. But in a suitable reference frame, which is an uni- 
formly rotating frame in the circular problem, and a non 
uniforml y rotating and pu lsating frame in the elliptic prob- 
lem (see ISzebehelvl l|l967l )). these periodic orbits become 
fixed points. When the motion of the secondary is circu- 
lar, these fixed points are elliptic (linearly stable) as long as 
the inequality 27m p mo < (m p + mo) 2 is satisfied l|Gascheaul 
1843). This ensures the linear stability of the equilateral 
equilibria as soon as h belongs to the interval [0, Ho], with 
Ho = (1 - \/23/27)/2 « 0.0385. For eccentric motions of 
the secondary, the (linear) stability criterion depends on its 
eccentri city, and stability exists for values of h greater than 
Ho f see (iDanbv 1964!; Roberts|[20o3 ') for numerical estimates 
and (|Mever fc Schmidtl2005h for analytical ones) . When the 
equilibrium is linearly stable, the eigenfrequencies (moduli 



of the eigenvalues) 



9l a 



and sl 45 , which yield the 



fundamental frequencies associated to libration, precession 
of the perihelion and precession of the ascending node of the 
asteroid, are given by: 



u Li5 = V(27/4)en p + o(y/i) = O(Vm), 
g L4:5 = (27/8)en p + o(e) = O(h) 



(1) 



withe 



m + m p 



0(h). 



General theory gives the values of the frequencies not 
only at the Lagrangian points, but also in the tadpole regio n 
as well as for horseshoe orbits (|Garfinkelll 19761 ; iMoraisl fcoO 1 ) . 
According to these theories, the magnitude of the frequen- 
cies remains unaltered in the whole phase space, except for 
s which reaches the order O(h) but remains always lower 
than g (in absolute value). Owing to these three different 
time scales, significant resonances between these frequencies 
are very unlikely (at least asymptotically). In the present 
Solar system, these resonances are negligible as long as h 
is lower than 10 -3 (close to Jupiter mass). Several regions 
of Jupiter's Trojan swarms are deeply affected by secondary 
resonances involving the commensurability of the libration 
frequency with Jupiter's proper mean motion. In Paper I, 
these resonances, elements of Family I , are defined by: 



iv + jn p + kg = 0, 



k) ez 3 



(2) 



For Jupiter's Trojan swarms, among this family, the in- 
fluence of the resonances satisfying j = 1 and i £ {12, 13, 14} 
is dominant for moderate to high libration amplitude. In the 
present paper, these three resonances can be seen in Fig. 
[5] (see section f2 . 3 [) . Although the dynamical implication of 



these resonances is weaker in the circular RTBP than in the 
elliptic one, Family I is well known in the first model. Indeed, 
the existence of denominators associated to the resonances 
€ {(11 : -1),(12 : -1), (13 : -1),(14 : -1)} dur- 
ing the Birkhoff norma lization process is already mentioned 
in (|Deprit et al.l Il967h . But these terms do not generate 
any difficulty up to degree 15. Moreo ver, those resonances 
are identified in (|Giorgilli et al.| [l989) as the factor of di- 
vergence o f the n ormal forms in the spatial circular RTBP. 
lErdi et all (|2007h study in detail these secondary resonances 
in the elliptic RTBP and their dependence on the mass and 
eccentricity of the secondary. 



2.3 Restricted (N+2)-body problem 

The first difficulty comes from the fact that, in the re- 
stricted (N+2)-body problem, since N is greater or equal 
to 2, the triangular equilibrium points do not exist. These 
points are replaced by elliptic qua s iperiodic traje ctories (see 
l|jorba fc Simdlll996l : [jorbal l2000l : iGabernl |2003| V) . But nu- 
merical experiments show that, if the additional planetary 
perturbations are not too large, a stable region remains in 
the neighborhood of the equilateral points. As in the RTBP, 
these equilateral points, which are always properly defined, 
are usually called L4 and L5. 

This point being clarified, let us study the influence of 
planetary perturbations on the Trojans associated to the p th 
planet. These dynamical effects on the Trojan swarms can 
be mainly split in two classes: the direct influences, due to 
the gravitational attraction of the planets on the Trojans; 
and the indirect effects, coming from the fact that the p th 
planet does not evolve any more on a Keplerian trajectorjjf]. 
Although this difference between direct and indirect pertur- 
bation is quite arbitrary (we will see later that in some cases 
these two phenomena are mixed together), it has the advan- 
tage to simplify the discussion. In terms of frequencies, the 
main effect of the direct influences is to modify the proper 
frequencies (y, g, s) of the Trojans. Assuming that the mo- 
tion of the additional planets is circular (we consider here 
the iV-circular problem), according to the Laplace-Lagrange 
theory, the secular linear contribution of the j th planet to 
the proper precession frequency g can be roughly approxi- 
mated at L4 (or L5) by the expression: 



a - - < I- 

a v 



(3) 



„ = in^a&'/ate). if a, = *± < 1. (4) 



9l, 



where b^L is a Laplace coefficient and Hi ~ rrij/mo (see 



for example dMurrav fc Dermott|[l999l )). Then, the secular 
frequencies at the triangular equilibrium points are given by: 



(Tot) _ (R.T.B) . U) 



J™) 



J.R.T.B) _ (j) 

I "^4,5 



(5) 
(6) 



2 In order to simplify the following discussion, we approximate 
all the motions as quasiperiodic, which is not necessarily the case. 
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where g^ L 



(R.T.B) 

4,5 



(27/8)en p and 



^4,5 



0, see formulas 



' ' and s' r To '' as well as g[ R ' T B ^ and 

4,5 a l>i.E 



{T} . The frequencies gr}^ 5 a™ o £ 

the planetary contribution = Sj^p 3l 4 ' 5 are gath- 

ered in Table ([T]l. for every planet of the Solar system. Al- 
though the results of Table ([TJ are valid only at the tri- 
angular points and for circular and coplanar planetary or- 
bits, they are sufficient to deduce striking conclusions. As we 
can see by comparison of the second and third column, the 
planetary contribution g^jV 1 ^ is not always a perturbation, 
but, at least for Mercury, Mars and Uranus, it can be the 
main contribution to the precession frequencies. It is only for 
Saturn and especially for Jupiter that the direct perturba- 
tions by the other planets impose only small corrections to 
their Trojans' secular frequencies. Consequently, except for 
Jupiter, the planetary direct gravitational attraction dras- 
tically modifies the dynamics predicted by the RTBP. This 
pheno menon has been emphasized by [Tabachnik fc Evansl 
(j2000|) in the case of Mercury's Trojans. Indeed, in the 
full Solar system model, the authors found that the most 
stable zones (100 Mys stability) did not contain the La- 
grange points. Another striking point is that, in an outer 
Solar system simulation, the region of smallest amplitude 
of libration of Uranian and Neptunian Trojans is shifted by 
more than 0.1 AU from it s predicted location in the RTBP 
ijNesvornv fc Donedliool ). 

The indirect perturbations act in a more subtle way 
than the direct ones: the introduction of forcing frequencies 
(i.e. the planetary frequencies, which are constant once the 
planetary system is given) enables resonances between the 
Trojan's frequencies and these additional frequencies. Obvi- 
ously, the secondary resonances with the orbital frequency 
n p are directly affected by the additional secular frequen- 
cies. These forcing frequencies increase the number of possi- 
ble resonant harmonics, imposing a generalization of Family 
I by: 



iu + jn p = —(kg + Is + kg 



Is 



(7) 



with i 7^ 0, j 7^ and j + k + 1 + k g ■ 1 + l s • 1 = 0. In these ex- 
pressions the dot denotes the Euclidian scalar product and 
1 = (1, • • • ,1). The enrichment of this family may generate 
large chaotic zones due to the overlapping of those reso- 
nances. Another expected consequence is the introduction 
of secular resonances. In the Solar system, as is shown by 
the last two columns of Table[TJ except for the frequency g of 
Jupiter's and Saturn's Trojans, the secular frequencies at L4 
are very close to the fundamental frequencies o f the perihelia 
and o f the nodes of the planets (|Laskarlll99(il : lLaskar et al.l 
2004). Indeed, the resonance s = S2 which limits the long- 
term stability regio n of low inclin ation Jovian Trojan^f] was 
first mentioned by lYoderl (11979m. other secular resonances 
involving s are discussed in (|Milanil Il993l Il994). The role 
of secular resonances in the moti on of the Trojans of th e 
inner planets was investigated by iBrasser fc Lehtol l|2002D . 
More accurate studies were dedicated to the role of secular 
resonances in the motion of Venus Tr ojans (Miche j 1 19971 : 
IScholl et al1l2005bl ) and Mars Trojans (jScholl et al.ll2005al ). 



3 Here, S2 has to be understood as the fundamental frequency 
associated to the precession of the ascending node of Saturn. This 
frequency is usually denoted s@. 



Table 1. Linear secular frequencies evaluated at the trian- 
gular points of the planets of the Solar system, expressed in 
arcsec yr -1 . In the first column is the name of the planets. The 
second one gives the precession frequency of the perihelion de- 
rived from the circular RTBP (formula HJ), while the third one 
gathers the contribution of the planets to this frequency (formula 
(|6]l). Finally, the precession frequencies g and s at L4 or L5 are 
given in fourth and fifth columns. 



Me 


(R.T.B.) 

3.0 


(plan) 

5.5 


(Tot) 

8.5 


s (Tot) 
-^4,5 

-5.5 


V 


17.3 


12.2 


29.5 


-12.2 


E 


13.3 


12.9 


26.2 


-12.9 


Ma 


0.7 


17.8 


18.6 


-17.8 


J 


352.3 


7.4 


359.8 


-7.4 


S 


42.3 


18.3 


60.6 


-18.3 


U 


2.3 


2.7 


5.0 


-2.7 


N 


1.4 


0.7 


2.0 


-0.7 



As in Paper I , the family containing secular resonances is 
denoted Family III and defined by: 



kg + Is + kg<Tg + l B cr s — 0. 



(8) 



Contrarily to Paper I , where k is always set equal to zero, 
the possibility that fc 7^ in Jupiter's Trojan swarms is 
discussed in sections !3.1.2l and l3.3l 

The conjunction of both direct perturbations and 
secular resonances prevents inner planets, except Mars 
(jScholl et al.l l200 5ai, from having long-lived Trojans. For 
this reason, the Trojans of the first three planets of the So- 
lar system are generally transient objects, which spend only 
a few hundreds of thousands yea rs in the co-orbital region 
IIMorais fc Morbidelli|[2002l . l2006l ). 

Until now, only secular forcing frequencies have been 
taken into consideration, but combinations of planetary 
mean-motions also play a major role in some specific con- 
figurations of the planetary system. If we consider a (5 : a 
MMR between the p th planet and another onjfl, let us say 
the q , its critical angle reads: 



= a\ p — p\ q + 



(9) 



where the dots represent a linear combination of longitudes 
of the nodes and of the perihelia such that the d'Alembert 
rules are satisfied. The modulus of the quantity v a ^ — 
an p — /3n q will depend on how close the planetary system 
is to the resonance. The closest the system will be to the 
resonance, the smallest \f a ,p\ will be. Hence, far from the 
MMR, v a! p — 0(1) (same order as the planetary mean mo- 
tions). The quantity \v a ,g\ becoming smaller and smaller 
as the two planets approach the exact resonance, this fre- 
quency will reach values that can generate resonances with 
the fundamental frequencies of the Trojans. As the plane- 
tary system approaches the MMR, |^ Q .^| will at first be of 
the same order of magnitude as v (i.e. 0(y / e)), enabling com- 
mensurabilities which generate the resonances of Family II , 
defined by: 

if - jv a ,p = - (kg + Is + kg ■ rr g + l s • <r s ) (10) 
with j(/3 — a) + k + / + kg ■ 1 + l s ■ 1 = 0. This generalizes 

4 MMR involving three and more bodies are not taken into ac- 
count here. 
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8.0 8.1 8.2 8.3 8.4 8.5 




9.50 9.52 9.54 9.56 9.58 9.60 9.62 9.64 



-9 -8 -7 -6 -5 -4 -3 -2 

Figure 1. Section of the phase space of the planetary system 
by the plane of initial conditions (a2,M2). Neighborhood of the 
2:1 (top) and 5:2 (bottom) mean motion resonances. The small 
bottom strip displays the grey code associated to the diffusion 
index: the darker the grey, the more regular the trajectories. The 
two horizontal white lines correspond respectively to the "elliptic 
segment" for M2 = 340.04° and to the "hyperbolic segment" for 
M 2 = 24.14° 

Family II as it is defined in Paper I . Once this threshold is 
crossed, no new significant resonance arises until the planets 
are very close to the MMR. Here, i/ a ^ is of the same order 
as g (i.e. 0(e)), which generates the resonances of Family 
IV satisfying: 

j^ Q ,/3 + kg = -(Is + kg ■ u g + l s ■ rr s ) (11) 

with j(a — f3) + k + I + k g ■ 1 + l s ■ 1 = 0. This generalizes 
the family Family IV that is denned in Paper I. 

A similar phenomenon arises after the crossing of the 
orbital resonance, when the frequency \v a ,/3\ increases from 
zero to O(l). 



3 APPLICATION TO JUPITER'S TROJANS 
3.1 One parameter model 

3.1.1 Model and method 0} analysis 

In this section, the evolution of the resonant structure in 
the Trojan swarms mentioned in section 12.31 will be illus- 



trated using a concrete planetary system. To this purpose, 
the methods described in Paper I will be followed. We con- 
sider the planetary system made of the Sun, Jupiter and 
Saturn. In order to be consistent with the previous section, 
the index 1 is associated to Jupiter and 2 to Saturn. In 
Paper I , the system, in its present configuration, is close to 
the 5:2 MMR. Consequently, resonances of Family IV should 
be present in the Trojans' phase space. Indeed, several res- 
onances of this family (for j = 1 and k — 4 in formula 
(|11[) ). have been clearly identified in Paper I . We have shown 
in this paper, that these resonances are associated to their 
long-term erosion of the Trojan swarm. 

Despite the long distance between the theoretical loca- 
tion of the 2:1 resonance and the present couple of giant 
planets (more than 1.2AU) resonance members of Family 
//associated to 1/2,1 have also been clearly identified. Owing 
to this distance, the combinations between v and fi,2 are of 
high order: i — 5 and j — —2 in formula (|10[) . Obviously, 
secondary resonances of Family / and secular resonances of 
Family III play a major role in Jupiter's present swarms of 
Trojans (Paper I). But, as these families are only weakly 
affected by migration (section 13. 1.2f> . except when they are 
very close to a MMR (section I3.3f> . we will not really pay 
attention to these resonances. 

In order to appreciate the modifications of the resonant 
structure of Jupiter's Trojans due to different relative posi- 
tions (and distances) between the two planets, we consider 
a sequence of independent planetary initial conditions. The 
fundamental frequencies of the planetary system depending 
mainly on the semi-major axes of the planets, or more pre- 
cisely on their ratio, we have decided to change only one 
parameter: the semi-major axis of Saturn 02. Hence, we are 
left with a one-parameter model. More concretely, for every 
value of the parameter 02, we integrate the Sun, Jupiter, Sat- 
urn and a set of fictitious Jupiter's Trojans (considered as 
test-particles). Except for the initial semi-major axis of Sat- 
urn, which is chosen between 8 and 9.7 AU, the initial con- 
ditions of Jupiter and Saturn are the ones given by DE405 
at the Julian date 245 2200.5 (2001 October 10). Regarding 
the Trojan swarms, the initial elements of their members are 
independent of the parameter. In each run, a grid of initial 
conditions is considered with 200 initial values of the semi- 
major axis and 40 values of the eccentricity equally spaced 
in the domain A = [5.2035, 5.4030] x [0.05, 0.30] (8000 test- 
particles) , while the other elements are fixed to the following 
values: a — A — Ai = 7r/3, <t 9 = zu — zu± — n/3, Q, = fli 
and I — Ii + 1* . Because the resonant structure depends on 
the initial inclination (Paper I), I* is fixed to three differ- 
ent values: 2° , 20° and 30° . In this paper, this set of initial 
conditions is denoted T>x* ■ 

The numerical simula tions are performed by the sym- 
plectic integrator SABA4 (|Laskar fc Robutell200ll ) with an 
integration step of 1/2 year. Trojans and planets are inte- 
grated on two consecutive time-spans of 5Myr. For particles 
surviving the integration (bodies which are not ejected from 
the co-orbital region before the end of the integration) the 
fundamental frequencies are computed for each of these two 
time-int ervals u s ing th e frequency analysis method devel- 
oped by lLaskaij (|l990h . If we denote by T the map which 
associates to ea ch Trojan of Z>z* its fundamental frequen- 
cies (u, g, s) (sec Laskar 1999 and Paper I ), the domain of the 
frequency space reached by the Trojans is: Oj* = J-(T>x*). 



6 P. Robutel and J. Bodossian 



Consequently, three types of complementary information are 
derived from the study of T>x* . in Paper I The most straight- 
forward piece of information is given by the escape rate: the 
number of Trojans escaping the co-orbital region before the 
end of the integration (10 My), divided by the initial num- 
ber of Trojans inside 2?i* . This indicator will be widely used 
in section [3T3T Owing to a loss of accuracy of our integrator 
during close encounters with Jupiter, the depletion of T>x* is 
probably over-estimated, but in any case, this ejection rate 
is always correlated to the global instability of the consid- 
ered region. Other significant information can be found in 
the frequencies. Since a detailed discussion of the applica- 
tion of Frequency Map Analysis to the Trojans can be found 
in Paper I , let us mention two applications of this method. 
First, the comparison of the set of frequencies computed on 
the two intervals of 5Myr makes it possible to derive the 
diffusion rate of every Trojan (index related to its stability). 
Practically, if we denote by i/ 1 ' the libration frequency de- 
termined on the first time span, and the same quantity 
computed on the second interval, the diffusion index will be 

(1) (2) 
V — V 

given by the quantity: log 10 | ^ . From this index, a 

dynamical map of the domain T>i* is derived. Second, the 
study of the frequency domain 0/* enables us to understand 
how the chaotic regions are generated by the overlapping of 
the underlying resonances. These complementary techniques 
will be used in section 13.2.21 

Once the fundamental planetary frequencies are known, 
it is straightforward to predict whether the families of reso- 
nances associated with these frequencies are inside the Tro- 
jans' phase space or not. Obviously, in order to make these 
predictions, the bounds of the frequency domain O/* have 
to be known. To this aim, we assume that these bounds do 
not depend on the value of 0,2, which is a very good approxi- 
mation in the case of Jupiter's Trojans. As it was established 
in Paper I , we assume in the following sections that: 

9 2 C [7700,9150] x [310,445] x [-45,-7.5] 

620 C [7400, 8660] x [285, 350] x [-40, -3.5] (12) 

630 C [7000, 8138] x [251, 280] x [-30, +0.6] 

where the three intervals are respectively the projections 
of &i* on the 1-dimensional space of u, g and s denoted 
n v (Qi*), 7r g (0/*) and n s (Oi*). The units of frequency are 
arcsec yr - . It is important to mention that the lower 
bounds of v and s here are arbitrary. For example, at 2° 
of initial inclination, the Trojans whose libration frequency 
is lower than about 7800 arcsec yr~ have trajectories rel- 
atively far from quasi-periodic, making the use of funda- 
mental frequencies less significant (details can be found in 
Paper I). The same remark holds when s is lower than the 
bounds indicated in (|12[l . These bounds may also b e com- 
pared to those obtained b y analytical fits of (v,g,s) (|Milanil 
1 19941 : iMarzari et alllioolh , which give quite similar results 
up to 30° of initial inclination. 

3.1.2 Behavior of the planetary frequencies 

From now on, Jupiter's and Saturn's semi-major axes and 
proper mean motions are denoted (01,02) and (711,712). Fol- 
lowing the notations established in section 12.11 the secular 
frequencies of the planetary system are cr g = (<?i,<?2) and 



cr 3 — (0, S2). These frequencies play a fundamental role in 
this study. It is therefore essential to know their variations 
with respect to the parameter 02. We have seen, in section 
12.31 that the location of the planetary system with respect 
to MMRs plays a major role in the transition from the reso- 
nances of Family II to the ones of Family IV . It is therefore 
important to know the geometry of the orbital resonances 
in the phase space. Fig. [1] shows dynamical maps of the re- 
gions surrounding the orbital resonance 2:1 (top) and 5:2 
(bottom), corresponding to the section of the phase space 
by the plane of coordinates (02, M2), M2 being the mean 
anomaly of Saturn. The grey code indicates the diffusion 
rate of Saturn's proper mean motion computed on two con- 
secutive time-intervals (section 13. 1 . 1 p . From these two de- 
terminations, denoted respectively and , a diffusion 
index is derived by the expression: log 10 \ (n^ — n^)/n^\. 
The structures of high diffusion rate (light grey for diffu- 
sion rate greater than —3) are associated to the inner part 
of the resonance (libration island) . Fig. Q] shows the char- 
acteristic shape of the resonant chains c omposed of two is- 
lands for the 2:1 and five for the 5:2 fsee lRobutel fc Laskarl 
(2001) for details). These large islands, which correspond 
to globally stable regions (elliptic regions in the pendulum 
model) , are separated by narrow unstable structures like the 
hyperbolic fixed points in the simple case of the pendulum. 
Even if the description of the dynamics of this problem is 
outside the scope of this paper, two points are interesting 
to note. First, the dynamics inside the resonant islands ap- 
pears very rich. Indeed, very sharp structures indicated by 
different diffusion rates are clearly visible; they are probably 
related to secondary or secular resonances. Second, particu- 
larly for the 5:2 MMR, the island chain is strongly distorted, 
making difficult to define a resonance width. The segment 
containing the initial conditions used for our simulation is 
represented in Fig. [T] Along this segment, the initial values 
of M2 are always equal to M2 = 340.04°. This line of initial 
conditions crosses the first lobe of the 2:1 orbital resonance 
nearly along its widest section, and passes near the libra- 
tion center. For this reason, we cal l it "elliptic segment". 
According to iMorbidelli et ail (|2005l ). during planetary mi- 
gration, the MMRs, and particularly the 2:1 are "jumped" 
by the planetary system. More precisely, if the migration is 
slow enough to satisfy the adiabatic invariance hypothesis, 
the system has to cross the resonance through its hyper- 
bolic fixed point (this makes sense for one degree of freedom 
systems) without reachi ng any libration zone (t he same phe- 
nomenon is observed in IMarzari fc Schollll2007h . In order to 
compare the "elliptic crossing" to the "hyperbolic" one, we 
consider a second set of initial conditions represented in Fig. 
□ by the second horizontal white line (M a = 24.14°). This 
path will be called "hyperbolic segment". We also notice 
in Fig[T] (bottom) that this segment passes through the hy- 
perbolic region connecting the lobes of the 5:2 resonance. 
Keeping this figure in mind, it is now easy to understand 
the behavior of the planetary frequencies under variations 
of <i2 along the two above mentioned paths. 

Jupiter's proper mean motion m is practically not af- 
fected by the variations of 02. As a result, the location of the 
resonances of Family I is almost independent of the value of 
d2- Only small shifts are observable due to the variation of 
the secular frequencies involved in formula 0. On the con- 
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Figure 2. Fundamental frequencies of the planetary system versus a<2- Top: combinations of frequencies 1/1,2!, 1^3, 7| an d 1^2 s| which 
are associated to the resonances of Family II . Bottom: secular frequencies gi,g2 and S2- The empty intervals correspond to the 2:1, 
7:3 and 5:2 MMRs where the secular frequencies are singular. The frequencies associated to the "elliptic" section of the phase space 
(M2 = 340.04°) are plotted in red, while the "hyperbolic" section [M2 = 24.14°) is in green. See the text for more details. 



trary, 72,2 being affected by the variation of the parameter, 
the mean motion combinations associated to orbital reso- 
nances (here 1/1,2, 1/3,7 and 1/2,5) vary drasticalljQ, and enable 
the resonances of Family II to go through the whole Trojan 
swarm. The evolution of the planetary frequencies along the 
elliptical segment is represented in Fig. [2] by red solid lines, 
while the green dashed curves show their evolution along 
the hyperbolic segment. Let us focus on the elliptic segment. 
Before going further, we have to mention that the quantity 
v a ,p is an increasing function of 0,2 as long as j3 > 0. Even 
if it is more convenient to plot \v a ,p\ rather than its signed 
value, we have to keep in mind that this frequency is nega- 
tive before an orbital resonance (i.e. when 02 is smaller than 
the value required to be in MMR with Jupiter) and positive 
after. When 0,2 evolves towards a MMR, the corresponding 
combination \v a ,p\ decreases towards zero. First, this fre- 
quency reaches values close to 8000 arcsec yr _1 giving rise 
to the resonances v w v a ,g of the second family. Then, as 
d2 increases, \v a ,p\ starts decreasing until it becomes com- 
parable to a few g, where some resonance of family IV is 
encountered. When the planetary system crosses the reso- 
nance, the frequency u a ,f3, which is equal or at least very 



Only these three resonances are taken into account in our 
study. Other MMRs of higher order like the 9:4, 11:5 and 12:5 
MMRs are also present in the domain of study, but their influence 
on the Trojans remains too small to be appreciated. 



close to zero, suffers from chaotic variations due to the dy- 
namical structures encountered inside the MMR (see Fig. 
[TJ. For the sake of clarity, the behavior of the frequencies 
inside orbital resonances is not reproduced in Fig. [5] This 
leads to gaps in the curves representing the frequencies. On 
both sides of these gaps, the singularities generated by the 
stable and unstable manifolds of the resonance impose that 
the frequencies go to infinity. 

The bottom frame of Fig. [2] displays the values of the 
secular frequencies (red curves) with respect to 02. The neg- 
ative frequency is S2, while the two positive frequencies are 
gi and 32 knowing that g\ < g%. The general trend of the 
absolute value of the secular frequencies is to decrease when 
a,2 increases. This is merely due to the fact that the per- 
turbations between the two planets decrease with respect to 
their mutual distance. But, as for the combinations of mean 
motions, singularities appear when the " separatrices" of the 
MMRs are reached. The increase in the secular frequencies 
is particularly striking on both sides of the 2:1 MMR. Here, 
as it can also be seen in a small neighborhood of the 7:3 and 
the 5:2 MMRs, <?2 is much more affected than the two other 
secular frequencies. More precisely, these three frequencies 
have to go to infinity, but the growth of the slope begins 
farther from the "separatrix" for 32 than for g\ and S2- As a 
result, we would expect to detect the influence of the secular 
resonances of Family ///like the g — (k + l)g2 — kgi for at 
least small values of \k\. But, as large values of g\ and g2 are 
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necessary to reach these resonances, this phenomenon oc- 
curs only very close to the 2:1 MMR. As it will be shown in 
section [3~3l instability is so strong in this region that secular 
resonances cannot be isolated and clearly identified. 

The behavior of the frequencies along the "hyperbolic 
segment" is qualitatively the same as along the "elliptic seg- 
ment", except near the three dominant MMRs, where the 
singularities are shifted as predicted by Fig. [T] Indeed, the 
approach of a MMR is characterized by a sharp and sud- 
den variation of the planetary fundamental frequencies cor- 
responding to the singularity associated to the separatrix. 
But the location of these asymptotes depends on the initial 
angles of the planets (M, oj,Q). The consequences of this 
dependence on the initial phases will be studied in section 
1331 



3.2 Sweeping of Family IV 's resonances across the 
Trojan region. 

3.2.1 Prediction of the location of Family IV 's 
resonances 

In this section we focus our attention on the close neighbor- 
hood of the 5:2 orbital resonance. Indeed, we have shown in 
Paper I that the system is close enough to this resonance to 
give rise to narrow unstable regions resulting from Family 
IV . Even if these chaotic regions seem very thin, several 
resonances of this family are involved in the low erosion 
pr ocess of Jupite r 's Tro jan swarms that was first mentioned 
bv lLevison et al.l (|l997T ). Because the Trojans will encounter 
several resonances of this family between the current loca- 
tion of the system and the 5:2 MMR, the study of this region 
enables us to describe and illustrate accurately the crossing 
of the Trojan swarms by the resonances of Family IV . 

The resonances of Family IV , defined by formula (111[) . 
involve the frequencies z/2,5, g, and the secular frequencies of 
the planetary system. According to Fig[3]and formula (I12p . 
the variations of gi,g2, S2 are small with respect to those of 
1/2,5. Furthermore, the secular frequencies of the planets are 
small compared to g, except very close to the 2:1 MMR (see 
later). As a result, the location of Family IV depends mainly 
on the values of ^2,5(02). The variation of this frequency 
with respect to 02 is plotted in figure [3J More precisely, 
the bold curve of this picture represents the graph of the 
function /: x 1 — ► — U2,s{x). The X-axis corresponds to the 
initial values of 02 (in AU) while the Y-axis is associated to 
the frequencies (in arcsec yr _1 ). This curve decreases very 
regularly, until 02 reaches 9.602AU, where a sharp change in 
the slope indicates the beginning of the 5:2 MMR (see also 
FigtU bottom-frame). It is now straightforward to predict 
the location of resonances belonging to Family IV in the 
Trojan phase space. Indeed, as long as the right side of pip 
is negligible with respect to kg, the resonant condition is 
well approximated by the relation: 

-^2,5(02)6^3(61.). (13) 
J 

For given values of j and fc, the previous formula defines a 
frequency interval where the resonance is reached. Approx- 
imating the bold curve in Fig. [3] as smooth and monotonic, 



/ _1 maps this interval in another one of the X-axifl namely: 
f~ l ((k/j)ir g (Gi*)). The cartesian product of these two in- 
tervals, the first one on the frequency axis (Y) and the sec- 
ond one on the 02 axis (X), defines a "resonant rectangle" 
displayed in gray in Fig|3] and denoted by R^'J . 

For the sake of clarity, our study is limited to the res- 
onances defined by j = 1 and k > (fc < on the other 
side of the MMR). Consequently, the rectangles R/* will 
be denoted Rj*. According to (|12l) and (|13[) . the location 
of the resonances of Family IV depends on the initial incli- 
nation 7*. For this reason, we study sections of the phase 
space at two initial inclinations: I* — 2° and 7* = 20°. The 
corresponding resonant rectangles are colored in light grey 
and dark grey respectively. The labels 5g, 4g, 3g and 2g in 
the bottom-left corner of the rectangles at 2° and in their 
top-right corner at 20° correspond to k = 5, 4, 3 and 2 re- 
spectively. Fig|3jshows that, for a given fc, the size of the Rj* 
decreases with 7*, and that for a fixed 7*, these rectangles 
are shifted rightwards as k decreases. These two phenom- 
ena are due to the facts that the function / is decreasing 
and concave, and that the width and the bounds of 7r ff (B/*) 
decrease while 7* increases. One of the main consequences 
of these properties lies in the fact that Family IV enters the 
Trojan swarms for a value of 02 which depends on 7*. There- 
fore, during a migration (slow monotonic variation of 02 with 
respect to the time in the studied region), Family IV's reso- 
nances will sweep through the Trojans phase space in a way 
that depends on the inclination: R% being larger than 7?2o> 
it follows that the instability generated by Family IV will 
destabilize the swarms much more efficiently at low inclina- 
tion than at higher inclination. Moreover, instabilities as- 
sociated to the fc-subfamilies of Family IV corresponding to 
Rj- , for different values of k, successively affect the swarms, 
a fact which further enhances the dependence of the degree 
of destabilization of the swarms on the value of the incli- 
nation 7* . Finally, as Fig|3j shows, the resonant rectangles 
can also overlap, ensuring the coexistence of at least two fc- 
subfamilies in the same swarm. Such intersections and their 
induced dynamics are studied in the next section. 



3.2.2 Crossing of the Trojan swarms by the resonances of 
Family IV . 

In order to verify these predictions and to illustrate the ge- 
ometry of the resonances of Family IV , we have integrated 
the 7/4 Trojan swarm for several initial values of Saturn's 
semi-major axis 02. The associated configurations of the 
planetary system are represented in Fig. [3] by vertical lines, 
and the initial values of 02 can be found in Tab. [2] The ini- 
tial inclination of the Trojans is equal to 2° for the solid 
lines, labeled with small roman letters, and to 20° for the 
broken lines labeled with capital letters. The present con- 
figuration of the Solar system, widely studied in Paper I , is 
not discussed here, but is represented in Fig(3] by the ver- 
tical bold line located at 02 ~ 9.5855^1(7. As ^2,5 is close 
to —1467 arcsec yr _1 , this implies that —^2,5 G 47r 9 (02) 



e Strictly speaking, the function / is not smooth on an interval, 
but, according to KAM theory, at most on a Cantor subset of this 
interval, which, however, is nearly of unit measure. 
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Figure 3. Crossing of the Trojan swarms by the resonances of Family IV . The slanted bold black curve represents the values of —^2,5 
( arcsec yr -1 ) versus 012 (AU). The grey rectangles indicate the values of (012, v) for which a subfamily of Family IV is inside the swarms. 
Light grey is used for I* = 2° and dark grey for 20°. The vertical bold line near 9.585AU corresponds to the present location of Saturn. 
The other lines indicate the values of ai used in the numerical simulations "a" to "f (I* =2°) and "A" to "F" (I* = 20°). See the text 
for details. 



(light rectangle) and that —^2,5 G 57r 9 (02o) (dark rectan- 
gle), which is in perfect agreement with the results of Pa- 
perl . 

The integer j is already fixed to 1, and from now on, we 
do not take into account the secular frequencies s and S2- 
This omission is fully justified for small inclinations, and we 
will see later that even for an inclination equal to 20°, the 
main resonances of Family IV are independent of s and S2. 
After this additional simplification, it is more convenient to 
rewrite relation pip as: 

V'l 5 k — 3 k 
9 = — + — — 9i + —(92 -91) (14) 

which satisfies relation when k = 3 — k gi — k 32 and 
k' = — k 32 . This new simplified formulation of the resonance 
condition defining Family IV enables us to exhibit some use- 
ful properties of these resonances. 

(i) For a given value of 112, the right-hand side of (|14|l is 
a constant denoted g r (for simplicity, the parameters k and 
k' are omitted). 

(ii) Still assuming fixed a^: as long as k and k' are not 
too big, the relation | 1*2,5 | 2> | (k — 3)gi + fc'(<?2 — gi)\ is satis- 
fied (see Section f3. 1. 2 [1 . Therefore, the resonance is reached 
for a value of g very close to — 1*2,5/ k. Consequently, we 
can see Family IV as being split in different subfamilies 
parametrized by k. In every subfamily, a single resonance is 
defined by the additional parameter k' . 

(iii) For a given k, the resonances of this subfamily are 
represented in the frequency space Oj* by parallel planes 
which are separated by a distance of ($2 — 9i)/k. Moreover, 
these planes are arranged in increasing order, in the sense 
that g r increases with k! . 

According to property (ii), the resonant rectangles 



provide an approximation of the location of the A;-subfamily 
which is valid as long as the combinations of the planetary 
secular frequencies involved in (|14p are small with respect to 
1*2,5/ k. Without going in too much detail, we can consider 
that these boxes are minimal in the sense that: at least one 
of the resonances of the fc-subfamily belongs to Rj* , and 
this resonance is always one whose dynamical influence is 
the largest. More precisely, except for k = 3, g = — 
does not satisfy equation (|14|l . and consequently does not 
correspond to any resonance of Family IV . But its closest 
resonant value, namely —V2,s/k + (k — 3)gi/k, is reached 
for k' = in formula (|14|) . If this resonant value does not 
belong to the box Rj, , one of the two frequencies — f2,s/k + 
(k — $)gi/k ± (#2 — gi)/k does. Furthermore, the resonance 
associated to k' = 0, or one of its two closest neighbors, 
plays a central dynamical role in the subfamily because its 
order is minimal (here we use the classical definition of order, 
that is: \k\ + \k — k' — 3| + \k'\), implying that its width are 
larger than the widths of all other resonances in the same 
family. Thus the condition g = — 1*2,5 /k suffices to predict 
approximately the value of 02 at which a major resonance 
of Family IV will influence the Trojans. 

Before going further, we will describe schematically the 
motion of a single resonance of Family IV during its cross- 
ing of the co-orbital region. Fig. [4] represents a section of 
the corresponding phase space in the (a, e) plane where 
the other initial elliptic elements are fixed. The horizontal 
straight line is the a-axis where e = e^ ' while the vertical 
grey line is the e-axis where a = af\ The two symmet- 
rical sides of the phase space section are represented here 
(a > and a < a^). Indeed, as it is explained in Pa- 
perl, the phase space of Jupiter's Trojans is symmetric. 
More precisely, we can clearly observe two different types 
of symmetries. The first one, which was already visible in 
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Figure 4. Schematic view of the displacement of a resonance of 
Family IV through the Trojan swarm. The X-axis represents the 
initial value of Trojan's semi-major axis. The initial eccentricity 
is associated to the Y-axis. The dashed curve approximates the 
limit of the stability region. The black curves correspond to the 
location of the same resonance for different values of Saturn's 
initial semi- major axis a.2- The larger the label, the larger ai- 

iMichtchenko et"al | (|200lr ); iNesvornv fc Donesl l|2002D . is a 
symmetry with respect to a curve that is close to the straight 
line a — ~ 5.2035 AU (i.e., the initial semi-major axis 
of Jupiter), and tangent to it at L4. The second symmetry 
is with respect to a curve close to the axis e = e^ ' ~ 0.0489. 
Moreover, considering quasiperiodic trajectories, the funda- 
mental frequencies corresponding to a given initial condi- 
tion and the ones corresponding to one of its two symmet- 
ric points are the same. These frequencies parametrize the 
KAM torus on which the given trajectories lie. This does 
not mean that the two corresponding trajectories are the 
same, but that they lie on the same invariant torus. From 
the dynamical point of view, these trajectories are equiv- 
alent. These symmetries point out the fact that there are 
manifolds (even clo se to L4) on which t he frequency map 
is degenerated (see iGabern et al.l l|2005t )), and they allow 
us to restrict the sample of initial conditions to the subset 
{(a, e); a > a^\e > }. 

The dashed grey curve gives the lower bounds of the 
strongly unstable region (the real form of this boundary cor- 
responds to the limits of the dark domains in all the right 
panels of Figs [5] and [6] which shall be discussed in detail 
below). Let us assume that we start the evolution of the 
planetary system with an initial value of 0,2 locating Sat- 
urn between Jupiter and the 5:2 MMR, in such a way that 
02 increases during the migration. As the planetary system 
gets closer to the 5:2 MMR, the frequency — 1^2,5 /k decreases 
towards zero ( 1/2,6 is negative in this region). Consequently 
the resonant frequency g r , associated to the fc-subfamily de- 
creases too, and the first contact between the Trojans and 
this resonance arises at the two black segments labeled with 
0, where a reaches its lowest and greatest value. Then, as 
g r keeps decreasing (and 0,2 keeps growing), the resonance 
goes towards the center (black curve 1) to reach L4 (label 2). 
Here, the two separated branches merge together to give a 
single curve. After the resonance travels through the vicin- 
ity of Z14, it moves towards higher eccentricities (label 3) 
and leaves the Trojan swarm through the secular resonance 
s = 82 (at least for low to moderate initial inclinations). 

This schematic view is essential to understand the evo- 



lution of the dynamics of Trojan swarms during planetary 
migration. Such evolution is presented in Figs[5]and|6] These 
figures show the sweeping of the co-orbital region, at I* = 2° 
and I* = 20° respectively, by the resonances of Family IV . 
Figs[5]and[6]are composed of two blocks. The right block cor- 
responds to dynamical maps of the domain D2 (Fig [5} and 
D20 (Fig [6]). As mentioned in section [3. 1.1 1 we use the rel- 
ative change of the frequency v given by a v = (i/ 1 — v 2 )/^ 1 , 
called diffusion index, as an indicator of the regularity of 
the motion. In this formula, v 1 is the libration frequency 
computed on the first 5 My, while v 2 is calculated on the 
following 5 My. Figs. [5] and [6] (right block) show dynami- 
cal maps in the action-like space (here a and e). A color 
is assign to each fictitious Trojan, coding its diffusion in- 
dex. The color scale ranges from blue, which corresponds to 
stable regions [a v < 10" 6 ), to red for very chaotic regions 
(ov > 10~ 2 ). In black, we display the particles that have 
been ejected from the Trojan swarms during the integra- 
tion (10 My). The left block is the corresponding view in 
the frequency space. The frequency map T establishes the 
correspondence between these two blocks which are dynam- 
ically equivalent. Consequently, figures [5] (left) are made of 
an union of curves (more or less smooth) which are the im- 
ages of the lines eo = constant by T . In the stable regions, 
the frequency map is very smooth implying the smoothness 
of the above mentioned curves (this is typically the case for 
s > —20 arcsec yr _1 in Fig. 0a and for s > —15 arcsec yr _1 
in Fig. [5]b). In contrast, in chaotic regions, T is singular, 
and the considered curves lose their smoothness. They are 
in fact disconnected and the points composing these curves 
seem to scatt er as in Fig. [5]a for s < —30 arcsec yr _1 (see 
lLaskarl (|l999h for more details). Each block is made of six 
panels, labeled from 'a' to 'F at I = 2° and A' to 'F' at 
/ = 20°, associated to different values of a 2 . These values 
are represented in Fig. [3] by vertical lines crossing the grey 
boxes. These values are chosen such that from 'a' to 'f ' (resp. 
'A' to 'F') Saturn's semi-major axis increases, allowing the 
fc-subfamilies of Family IV to cross the Trojan swarms fol- 
lowing the rules established in section 13, 2. II The resonances 
of Family /Fare clearly visible in the frequency maps (left 
blocks) as vertical structures made of gaps or of accumula- 
tions of dots. These are enhanced by vertical dashed lines 
for fc' 7^ and by solid lines for k' = 0, where k! is the 
integer which appears in formula (j 14|) as the parameter of 
the fc-subfamily. On the action side (right blocks), it is not 
so easy to identify these resonances. Although having the 
typical shape drawn in Fig|4] confusions are always possi- 
ble since the elements of Family II also have a quite similar 
form. The main resonances of Family IV , which are drawn 
in Figs [5] and [6] are gathered together in Tab. [2] The initial 
values of Saturn's semi-major axis are listed in its second 
column. The corresponding labels are in the first column. 
The two last columns correspond to the values of the inte- 
gers k and k' appearing in formula (114[l and which define 
the resonances of Family IV associated to the 2:5 MMR. 

The secular resonance s = S2 is also easily identifiable 
in the left block by a horizontal line corresponding to an 
accumulation of dots at s = — 26 arcsec yr -1 surrounded by 
gaps. This resonance splits the frequency space in two differ- 
ent regions: the part located above the resonance, which is 
globally regular assuming that the planetary system is not 
too close to the 5:2 MMR (it is true at least between the 
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Figure 5. Destabilization of the Trojan swarms by the resonances of Family IV at I* = 2°. Left block: Dynamical maps in the frequency 
space (projection on the (g,s) plane). The vertical lines show the predicted location of the resonances of Family TV involved in the 
dynamical process (see Table [2); while the horizontal line emphasizes the secular resonance s = sg. Right block: Dynamical maps in the 
action space (the initial conditions of the fictitious Trojans are chosen in the (a, e) plane, the four other elliptic elements being fixed). 

The diffusion rate defined by log 10 | - ^j- - | is coded by colors that vary from —8 (stable) to —2 (strongly chaotic). Each block is split 
in six panels labeled from "a" to "f corresponding to the numerical simulations described in section r3.2.2l See the text for more details. 
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Figure 6. Same as Fig[5]for I* = 20°. 
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Table 2. Resonances of Family /V represented in Figs [5] and 
[6] The two first columns indicate the label of the corresponding 
panel and the initial value of d2 used for the associated numerical 
simulation. The two last columns give the values of the integers 
k and k' defining jn formula 1 1 41 ) the resonances represented by 
vertical lines on each panel of Figs \5\ and \E\ 



Label 


a 2 (AU) 


k 


k' 




a 


9 


58986 


i 


-I,... 


3 








3 


-2, ■ • ■ 


4 


b 


9 


59364 


4 


2,3 










3 


-2, ••• 


4 


c 


9 


59513 


3 


-2, ••• 


4 




9 


59513 


2 


-5 




A 


9 


59630 


4 


2,3 






9 


59630 


3 


-3, •■• , 


2 


d 


9 


59684 


3 


-2, ■ ■ ■ , 


4 








2 


-5, ■ • ■ , 


-2 


B 


9 


59756 


3 


-3,--- 


3 


C 


9 


597816 


3 


-3, ••■ 


3 


D 


9 


59882 


3 


-3, ••■ 


3 


c 


9 


60004 


2 


-5, ■ • ■ 


2 


f 


9 


60057 


2 


-5, ■ • ■ 


3 


E 


9 


60088 


2 


-6, ■ • ■ , 


-2 


F 


9 


60197 


2 


-5, ■ • ■ 






Saturn's present location = 9.5855^1)7 and 9.595AC7 sim- 
ulation c), and the part lying below this resonance which is 
strongly chaotic. This secular resonance appears in the ac- 
tions space (right block) as the red arch crossing the vertical 
axis at about e = 0.25 and the horizontal one at a = 5.35AU . 
As mentioned above, Fig. [5]a reveals that the domain lo- 
cated inside this arch (blue to green color) is much more 
regular than the outside part (orange to black). For small 
eccentricities, s — S2 overlaps with the Family /resonance 
defined by formula ([7]) with i = 1 et j = 13 (red V-shape 
around a = 5.35AU) , generating this large chaotic region 
around 5.34AU. Family I is also acting around 5.4AU where 
the narrow red strip corresponds to i — 1 and j = 14. But, as 
in this experiment the semi-major axis of Saturn varies only 
very little (about one hundredth of AU (see Tab. [2}) only 
the resonances of Family IV move across the Trojan swarm, 
while those belonging to the three other families are practi- 
cally fixed. This makes it easier to recognize the resonances 
of Family IV . 

According to Paper I , the present system (vertical bold 
line in Fig. [3J is such that the 4-subfamily of Family IV is 
acting on low inclined Trojans. For slightly larger 02, both 
the 4-subfamily and the 3-subfamily are involved in the 
Trojans' dynamics. In Fig[5]a, the 4-subfamily is exiting 
the phase space (the corresponding resonances are located 
around g — 330 arcsec yr" 1 ), while the 3-subfamily pen- 
etrates inside. For k — 4, the subfamily mainly appears 
above the arch s = S2 at high eccentricities (e > 0.25), 
and strengthens the instability induced by the secular reso- 
nance. The 3-subfamily appears at the opposite side of the 
frequency space and is very difficult to identify in the ac- 
tion space because it is located in a region already filled by 
the resonances of Family /mentioned above (red V-shape 
surrounding 02 = 5.35 AU). It is worth mentioning that the 
narrow gap at a = 5.38AU in Fig. [S]a is generated by the 
Family IV resonance defined by k = 3 and k' = 1. This res- 
onance will play a major role in the next experiment. When 



C12 is slightly larger, as in simulation 'b' ( 02 = 9.59364AI7), 
the situation is much more interesting. Indeed, the subfamily 
k — 3 begins to reach the central region, generating strong 
instability associated to a rapid escape of numerous Tro- 
jans, particularly for a > 5.31AU. Another important point 
lies in the fact that this subfamily encounters elements of 
Family II . This overlap generates the deep gap located at 
5.31AU. The involved resonances of Family II are defined by 
the multiplet (a, f3, i,j, k, l,l s ) = (1, 2, 5, 2, 0, 0, 0) in formula 
(|10|) and are identifiable in Fig. [5]a by the nearly vertical 
orange structure between 5.31 and 5.32AU (Its location is 
practically unchanged from Fig. [S]a to Fig. [S]f ) . 

In panel 'c', the 3-subfamily continues to cross the Tro- 
jan swarm, and now occupies its core. The resonance asso- 
ciated to k' = is now very close to L4 (curve 2 in Fig 
[4j while previous ones (k' < 0) are on the other side of 
Li (curve 3 in Fig [4} . The overlap of these resonances (for 
k' 7^ 0) generates the huge red region above the resonance 
k' = 0, while the other part of the phase space begins to 
be a bit more regular. In panel 'd', the 3-subfamily moves 
outwards but its resonances are still acting for e > 0.13. 
The sharp transition between the black region and the blue 
one is due to the resonance defined by k = 3 and k' = 2. 
Below this resonance, the size of the stable region (blue to 
green) is increasing again. Let us remark that the effect of 
the 2-subfamily (right side of the frequency map) is not yet 
noticeable. Finally, when the 2-subfamily penetrates deep 
inside the Trojan swarm, the dynamics becomes much more 
chaotic. In panel 'e', a huge number of Trojans escapes the 
co-orbital region in less than ten million years, while in 'f, 
only small islands of temporary stability remain. The com- 
parison between 'f '-left and 'f '-right shows that each island is 
separated by two consecutive resonances of the 2-subfamily. 
This is possible because the distance between two resonances 
of the previous kind is larger than the one separating the res- 
onances of the 3 and 4-subfamilies. The simulation 'f ' is the 
last one of the sequence. Indeed, if Saturn get closer to the 
5:2 MMR, the entire Trojan swarm is rapidly cleared. At 
9.602AU, (see Fig.[3J), the planetary system is getting closer 
to the MMR, and the frequency —1*2,5 decreases abruptly. 
During this decrease, the 1-subfamily is reached, but at the 
same time, the planetary system crosses the 'separatrix' of 
the 5:2 MMR. The chaos induced by the separatrix crossing 
combined with the resonances of Family IV make the whole 
Trojan swarm strongly unstable. We will see in section 13.31 
that this phenomenon stops when the planetary system is 
inside the stable domain of the libration area associated to 
the 5:2 MMR. 

At an initial inclination of 20°, the mechanism is prac- 
tically the same, but a bit delayed. Indeed, in simulation 
'c' (/* = 2°) and 'C (/* = 20°) the dynamical situation is 
basically the same, but 'c' occurs at 02 = 9.59513^4)7 and 
'C at 9.597816AC7 . This is a direct consequence of the fact 
that the upper bound of 7r 9 (Or*) is a decreasing function 
of /*. This may have implications on the final distribution 
of the Trojans versus inclination, as we will see in the end 
of the present section. Looking at Fig. [5] (left block), one 
realizes that the main chaotic structures are vertical. This 
confirms the choice, which could seem arbitrary, that we 
have made in formula (|14p . Indeed, the resonances of Family 
IV which generate the most chaotic behavior do not depend 
on s,Sj. Conversely, these secular nodal frequencies are in- 



Resonant structure of Jupiter's Trojans-II 13 



volved in secular resonances. In addition to the resonance 
s — S2, which is still visible at 20°, two other resonances 
can be seen in Fig. [6] (right block). These structures form 
the two yellow thin arches (particularly in 'A') cutting the 
X-axis at 5.26 AU for the 3s — S2 — 2gi = (inner arch), 
and 5.29 AU for the 2s — 3gi + g% = (outer arch). As in 
Fig. [5] simulations 'A' to 'F' show the shifting of the reso- 
nances of Family IV . In 'A', the 4-subfamily is exiting the 
phase space while the 3-subfamily has already reached its 
center. The two resonances in both sides of L4 are defined 
by k — 3 and k' — —3, —2. They are associated to the two 
quasi-straight lines displayed in yellow-orange colors, which 
cross the two above-mentioned resonant arches. In 'B', the 
resonances are shifted leftwards (frequency space) and their 
overlap generates global chaos above e = 0.12 (location of 
the resonance k = 3, k' = — 1). Then, the resonance k' = 
reaches L4 (in 'C') and generates the sharp transition be- 
tween the blue region and the red region, while the size of the 
regular region keeps shrinking. Next, as the 3-family moves 
outwards, the stable area begins to increase up to simulation 
'E' where resonances of the 2-subfamily with large \k\ pen- 
etrate the phase space to finally eject most of the Trojans 
in 'F'. As in simulations performed with an initial inclina- 
tion of 2° , only small islands of temporary stability remain 
and are surrounded by gaps generated by the resonances of 
the 2-subfamily. The increase in the size of the stable re- 
gion between the crossing of two consecutive subfamilies of 
Family IV is something that we have not observed at low 
inclinations. The explanation of this difference is very sim- 
ple. We know from section 13.2. II that the rectangles R 2 o are 
larger tan the -R20 an d that the boxes R 2 o overlap while 
the R200 do not. It turns out that in the gap lying between 
J?2° an d ^20° i t ne Trojan swarm recovers, at least partially, 
its stability. This phenomenon can be also observed in Fig. 
[8jb around a — 9.6AU, where the fraction of ejected Trojans 
decreases suddenly from 0.8 to less than 0.6 to finally jump 
to 0.9 and more (see section [3731 for more details). 

To complete this section, let us point out that the above- 
mentioned evolution may lead to a lack of weakly inclined 
Trojans. Indeed, if we assume that the migration stopped 
close enough to the MMR, the sweeping of the Trojans would 
have taken place for small inclinations while it wouldn't have 
begun for higher inclinations. Unfortunately, in the Solar 
system, Jupiter and Saturn are too far from the 5:2 MMR 
to have generated this phenomenon; except if this system 
was, in the past, closer to the resonance than it is now. But 
it does not seem to be a realistic scenario. 
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Figure 7. The 3:7 MMR and its Family IV : ©2 projected on the 
(g, s) plan for £12 = 9.1833AC7. From right to left the six verticals 
lines correspond to the resonances 2g = —^3,7 — 2gi + k (92 — 9i) 
for — k £ {1, 2, 3, 4, 5, 6}. The frequency unit is the arcsec yr~ 1 . 



Table 3. Secular resonances: values of 02 (in AU) for which 
the three first resonances of the form g = 92 + k(g 2 — gi) en- 
ter the intervals O2, 020 , and ©30 for the "elliptic section" 
(M 2 = 340.04°). 



2° 
20° 
30° 



a = 92 

< 8.39 

< 8.395 
< 8.4 



292-51 g = 3g 2 - 2gx 



[8.41,8.425] 
[8.425,8.43] 
[8.435,8.44] 



[8.43, 8.455] 
[8.45,8.46] 
[8.47, 8.475] 



Trojans frequency domain O2 on the plane of coordinates 
(g,s). For this planetary configuration, which is very close 
to the 7:3 MMR, the combination of proper mean motion 
^3,7 = 3m — 7n2 is close to —423 arcsec yr _ , implying the 
relation —1/3,7/2 £ n g (Q 2 ). Consequently, several resonances 
of the 2-subfamily, defined by the relation: 



2g = -1/3,7 -2gi +k 



•91) 



(15) 



are present in the Trojans' phase space. Fig. [7] shows those 
resonances for k' £ {— 6, — 5, — 4, — 3, — 2, — 1}. Let us notice 
that the horizontal line, corresponding to the secular reso- 
nance s = S2, is shifted downward in comparison with Fig. 
[5] Indeed, S2 was equal to about —26 arcsec yr _1 in the 
simulations of section [3.2.21 while here, its value is close to 
—32 arcsec yr _1 . 



3.2.3 Family IV associated to other orbital resonances 



Since the mechanism studied in section l3.2.2l is very general, 
resonances of Family IV should take place in a small neigh- 
borhood of each MMR. But in some cases, the corresponding 
region is so chaotic that no resonance can be accurately iden- 
tified. This is particularly true for the 2:1 orbital resonance, 
where its merging with Family II generates a huge unstable 
zone (see section EPfl . 

On the other hand, the resonances of Family 
IV associated to the 7:3 MMR can be easily detected, es- 
sentially because they generate moderate chaos. We con- 
sider here that the initial semi-major axis of Saturn is equal 
to ci2 = 9.1833AI7. Fig. [7] displays the projection of the 



3.3 Role of the secondary resonances of Family II 

Up to now, we focussed on the resonances of Family IV to 
clearly illustrate the displacement of instabilities inside the 
Trojan swarms during planetary migration. But investi- 
gations regarding this family, which only acts on a rela- 
tively small neighborhood of its associated MMR, confine 
our study to specific periods when the planets are close to 
MMRs. In this section we will adopt a more global point of 
view by trying to find out, in the most exhaustive possible 
way, which are the resonances able to affect Trojan swarms 
during a migration phase. As mentioned in section 13.1.1 1 we 
limit the migration path to the domain a 2 £ [8, 9.65]A(7, 
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Figure 8. Chaos generated by the resonances of Family II: on the bottom frame (d), the frequencies fc'i,2,l / 3,7,f2,5 and some of their 
rational multiples are plotted versus a 2 . The rational numbers written on the curves (ij 'j)vt at g) correspond to the value of i/j, the values 
of (a, (3) being equal to (1, 2) for the red curves, (3, 7) for the black bold ones and (2, 5) for the green dashed curves. The horizontal 
grey, pink and blue strips indicate the reachable values of v in ©2, ©20 > an d ©30 • The three upper plots represent the relative number 
of ejected trojans n e j ec t (escape rate) during our numerical integration (section 13.1. 1> of respectively £>2° (panel c), T>20° (b) and X>3o° 
(a) for a2 £ [8.3,9.65] AU. The arrows indicate the location of the secular resonance g = 3g 2 — 2pi. 



which locates Saturn between the 2:1 and the 5:2 MMRs. 
Our goal can be reached in, at least, two different ways. The 
first one consists in the prediction of these events knowing 
on one hand the behavior of the planetary frequencies, and 
on the other hand the four families of resonances generating 
chaotic behaviors. The second way relies on numerical simu- 
lations of the same kind as the ones presented in section [3. 2. 2 1 
and 13.2.31 We have followed independently these two differ- 
ent paths and have found results in very good agreement. 
These results are presented in Fig [5] Among the numerous 
data produced by these numerical simulations, we only plot 
the ejection rate, that is: the number of escaping Trojans 
divided by the amount of bodies present in the initial popu- 
lation. Fig. [8] is split in four panels. The bottom frame ([8]d) 
shows the predicted location of the resonances generating 
instabilities, described in detail later on. The three other 
panels display the depletion rate (Y-axis) with respect to 
the initial semi-major axis of Saturn (X-axis, in AU). They 
correspond to simulations with initial inclination equal to 
/* = 2° in Fig[8]c , I* = 20° |[8jb), and I* = 30° ©a). 
On these depletion curves, the succession of steep peaks of 
ejection (local maximum of ejection) and deep troughs of 
stability (local minimum of ejection) provide useful indica- 



tions on the global dynamical behavior of the Trojan swarms 
along the migration path. 

Before going further, let us notice that, in the present 
Solar System, the Jovian Trojans are located in a trough of 
strong stability (02 = 9.5855AU). It may be amazing that 
the migration would stop where the stability of the Trojan 
swarms is maximal. But in fact, other planetary configura- 
tions leading to stable swarms exist. According to Fig.[S]a-c, 
the most stable regions lie between 9A3AU and the left edge 
of the 5:2 MMR. In particular, a 2 = 9.4T4C/ provides a nice 
solution to this problem which is perhaps the most stable 
that we have explored, particularly at an initial inclination 
of 20° and 30°. 

The ejection peaks present in Fig. [8] are obviously asso- 
ciated to resonances. We will show in this section that, except 
three peaks, the other ones are generated by resonances of 
Family II . Let us first consider these three particular cases. 
Those unstable regions are indirectly generated by the three 
main MMRs present in the studied domain: the 2:1 MMR 
for a 2 < 8AAU (in fact this is more a plateau than a peak, 
in this whole region the depletion rate is about 100%), the 
7:3 MMR at a 2 = 9.2AU and the 5:2 MMR for a 2 around 
9.6AU and more. As mentioned in section 13.1.2 1 both sides of 
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these peaks are generated by action of Family IV conjugated 
to the chaos introduced by the planets crossing the separa- 
trices of the MMRs. When the planetary system is deep 
inside the 7:3 and more importantly 5:2 MMRs, the deple- 
tion rate decreases suddenly, which suggests the existence 
of regions harboring stable Trojans. This is especially strik- 
ing at I* = 2°. This ph enomenon has been pointed out in 
iMarzari fc Scholll (|2007l ) where authors find long-lived Tro- 
jans when Jupiter and Saturn are in 2:1 orbital resonance. 
We do not observe the same stability in this resonance, as 
we can see in Fig(8] where the ejection rate is always close 
to 100%. These different behaviors can be explained quite 
easily. Indeed, Trojan stability depends strongly on the lo- 
cation of the planetary system with respect to the different 
topological structures of the considered MMR. If the plan- 
ets inside the resonance are close a separatrix (hyperbolic 
manifold) of one of these structures, the chaos induced by 
the planets may destabilize the whole Trojan swarm. On the 
contrary, if the planets lie in a stable region corresponding 
for example to an elliptic libration centre, their motion is in 
general close to quasiperiodic. Consequently, the Trojans do 
not suffer from chaos induced by the planetary system. But 
the motion of the Trojans is not necessarily regular. Sec- 
ondary resonances between v and the libration frequency 
associated to the 1:1 MMR can arise, leading to the ejection 
of numerous Trojans. Without studying the MMRs one by 
one, we cannot know which case corresponds to a given sit- 
uation. We can only mention that, around 9.63AU, where a 
region of relative stability takes place inside the 5:2 MMR, 
the angle 2Ai— 5A2+3ci72 is librating, and according to Fig.[U 
our 'elliptic migration path' crosses some kind of libration 
center (one of the most stable regions encountered in the 
section of the phase space). On the contrary, when Jupiter 
and Saturn are in 2:1 MMR (a 2 < 8AAU), almost all Tro- 
jans are ejected from the co-orbital zone. This is because 
our migration path always stays quite far from the libra- 
tion center of the 2:1 orbital resonance (see Fig. [T]). Fig(9] 
shows a very different behavior when the segment of initial 
conditions crosses the libration center. This figure shows the 
ejection rate along a path where the initial value of M2 is 
always equal to 316° (solid line with white circles) while the 
"elliptic segment" (A/2 = 340°) is represented with dashed 
lines with black circles. While the ejection rate along the 
elliptic segment is almost maximal, the solid line indicates 
the presence of a more stable region. Indeed, an increase of 
stability, denoted by a sharp change in the stability index 
(with a minimum of about 0.6) arises when 02 belongs to 
the interval [8.28 : 8.33] AU, which fits very well with one of 
the most stable regions (some kind of libration center) inside 
the 2:1 planetary MMR (see Fig. Qjj. It turns out that the 
possibility to find stable Trojans when the two planets are in 
MMR depends strongly on the geometry of this resonance. 

In order to demonstrate that, except the three above- 
mentioned peaks, the main unstable structures (ejection 
peaks) are generated by Family II , we will develop three 
different arguments. First, a straightforward reasoning will 
show that the three other families cannot induce these in- 
stabilities. The second argument will be based on the very 
good agreement between the simulations presented in Fig. 
[S]a-c and the predictions regarding the location of the res- 
onances of Family II . The last argument will consist on a 
detailed study of the region lying between 9.26 and 9AAU. 
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Figure 9. Two different crossings of the 2:1 MMR. The dashed 
line with black circles corresponds to the " elliptic section" already 
presented in Fig. [8] Along this line, the diffusion rate (vertical 
axis) is always close to 1. On the contrary, the solid line which 
crosses the 2:1 MMR along the line defined by M2 = 316° indi- 
cates the presence of most stable Trojan swarms for a 2 € [8.28 : 
8.33] AU (see the text for more details). 



First, we know that Family TV acts only on a small 
neighborhood of a given MMR. This prevents this fam- 
ily from producing instabilities like the one generating the 
largest and highest peak of Fig. [8]c. Regarding Family I , its 
resonances depend mainly on n\ which varies only very little 
in our simulations. Consequently, the resonances of Family 
/have approximately the same strength along the migra- 
tion path and thus cannot be responsible of these successive 
stability changes. The same argument holds for the secu- 
lar resonances of Family III . The secular resonance most 
affecting Jovian Trojans is the s — S2 resonance. It moves 
continuously (except in a small neighborhood of the MMRs) 
without coming close to L4. Other higher order secular res- 
onances involving nodal secular frequencies can eventually 
move inside the Trojans phase space, but generate only 
a very local chaos (see section 13 . 2 . 2 p . A very interesting 
point is connected to the existence of the secular resonances 
g — (k + l)gi — fc(?2, which can play an important role in in- 
creasing the eccentricity of the Trojans belonging to. On the 
base of these values of gi and 32 given in Fig.[2]and according 
to formulas (|12p . the exact location of these resonances can 
be derived. Tab. Ogives the locations of the main resonances 
of this kind. According to Fig. [8] these values correspond to 
a region which is close to the 2:1 MMR, and where very 
high peaks are present. We come back to this point later, 
but in any case, the influence of secular resonances seems to 
be very local and cannot generate this succession of highly 
unstable regions. 

Having excluded the significance of other families for the 
peaks of the ejection rate in Fig. [8] only Family //remains. 
This is not very surprising for at least two reasons. The first 
one is that we have shown in section 12.31 that these reso- 
nances were able to act far away from their associated MMR. 
In fact, Paper I reports that, in the present Solar system, nar- 
row unstable regions resulting from Family II appear in Qi* . 
The second reason lies in the fact that in iMorbidelli et alj 
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l|2005l ). the authors point out two strong instabilities which 
enable the capture of the future Jupiter Trojans. Accord- 
ing to the authors, these instabilities correspond to regions 
where v ~ 3i/i,2 and v ~ 2v\^- Obviously, these resonances 
are members of our Family II . In order to make sure that 
these peaks are really due to Family II , we will compare, 
as in section 13,2.21 our predictions to the simulations for 
which results are plotted in Fig. [8] To predict the loca- 
tions of Saturn where resonances of Family II penetrate the 
Trojan swarms, we proceed as for Family IV: neglecting 
the right-hand side of equation (|10[) . that is, the combina- 
tion of secular frequencies, the resonance condition reads: 
v = jv a ,p- The question is now to know whether or not, 
for a given <J2, the frequency \v a: p belongs to 7r„(Oj-»). The 
bottom frame of Fig. [5] answers the question: for the three 
main MMRs, namely (a, 13) G {(1, 2), (3, 7), (2, 5)}, the val- 
ues ^v a ,0 are plotted versus a^. The X-axis represents a 2 
in AU while the vertical axis measures the frequencies in 
arcsec yr _1 . A color is associated to each MMR, the curves 
^v a ^ (which are practically straight lines) are drawn in red 
for the Family //connected to the 2:1 MMR, in black for 
the one associated to the 7:3 MMR and in green for the 5:2 
MMR. The black labels on these lines indicate the values of 
j associated to the corresponding resonance. For the sake of 
clarity, these labels are aligned. The upper line (located at 
about 8600 arcsec yr _1 ) corresponds to 1/1,2, the second line 
(at 8200) to 1/3,7, and the last one (at 7600) to 1/2,5. The col- 
ors of the background correspond to the values of the initial 
inclination /*. Grey is associated to /* — 2°: more precisely, 
in the grey region, the frequencies (Y-axis) range in the in- 
terval ^(02°)- In the same way, in the pink rectangle, the 
frequencies belong to 7r„ (820° ) , and to 7r„ (630° ) for the cyan 
domain. The use of these different colors allows to predict 
easily the location of resonances for the three inclinations 
used in this paper. 

As was the case for Family IV , the criterion j^a,^ £ 
7Ti/(0/*) does not correspond to a single resonance in the 
Trojans phase space but rather to a multiplet of resonances 
associated to the linear combination of secular frequencies 
arising in the right hand side of (| 10p . Consequently narrow 
strips should replace the lines. But for the sake of simplicity, 
we prefer to keep the lines. 

In FigHa-c, between the 2:1 MMR (left) and 8.8 AU, 
three high peaks stand out for all inclinations (the ejection 
rate is practically equal to 100%) . As mentioned above, the 
resonances of Family //associated to the 2:1 MMR dominate 
this region; therefore, these three strong depletion zones are 
generated by these resonances. Indeed, these peaks coincide 
perfectly with the lines labeled with 3, 2 and 1 (bottom 
frame) associated to the relations: v equals 3, 2 and 1 times 
1/1.2 ■ These maxima of ejection match very well the predic- 
tion: the bases of the peaks coincide with the projection 
on the X-axis of the intersection of the corresponding line 
with the colored region (i.e. pink for /* = 2°). Moreover, 
the peaks are shifted rightwards when the inclination grows, 
which is due to the fact that red lines have positive slopes. 
We have to mention that it is not the first time that thes e 
structures are observed. Indeed, in iMorbidelli et al.l <|2005h . 
the authors have widely studied the region corresponding 
more or less to the interval [8.3, 8.65] AU in our simulations 
(see next section). They found two regions of strong deple- 
tion corresponding to our first two peaks. Moreover, in Fig. 



1 of the above-mentioned paper, the curve indicating the 
fraction of population that survives for 2 10 s My in the co- 
orbital region possesses a small singularity (just after 1.5 My 
of integration) which is not mentioned in the paper. Looking 
at Fig. [8] we clearly see that this singularity is generated by 
the subfamily v = §1/1,2- 

The main peak corresponding to the encounter of the 
se condary resonanc e v = 1 /1,2, is outside of the region studied 
in IMorbidelli et all (|2005h . but it seems that an important 
phenomenon occurs here. Indeed, the rightwards shift of this 
structure when /* increases is more pronounced than in the 
two previous peaks. This happens essentially because the 
slope of the red curve representing the frequency v = 1/1,2 is 
smaller than the ones of th e curv es v = 2i/i,2 and v — 3i/i,2. 
F. Marzari and H. Scholl (120071 ) mention the existence of 
this resonance but observe that its influence is weaker than 
that of the previous secondary resonances (the 1:2 and 1:3). 
It is not what we get in our study, but the two simula- 
tions are not easy to compare. The first simulation takes 
into account a forced migration, imposing Trojans to cross 
a given resonance quite rapidly. While in our case, the mi- 
gration being frozen, the resonances have enough time (here 
10 My) to eject the test-particles. In any case, we probably 
overestimate the influence of the resonances, but our goal is 
not to have a realistic simulation of Jupiter's Trojans dur- 
ing planetary migration, but to point out the regions or the 
events that are relevant in term of ejection (or injection) of 
bodies in the co-orbital region. In the same paper Marzari 
and Scholl stress the dominant role of the secular resonance 
g — (72- In our simulation, the region where g = g2 occurs is 
so close to the 2:1 MMR (see Table [3} and so chaotic that 
the fundamental frequencies are meaningless: Trojan trajec- 
tories are too far from quasiperiodic ones. But, on the other 
hand, for /* = 2° , we can distinguish a very small peak cen- 
tered at Gt2 = 8A4AU (arrows in Fig(8]a-c) which is gener- 
ated by the secular resonance g = 3(?2— 2<?i. By increasing /*, 
the secular resonance is shifted rightwards (Table [3} while 
the secondary resonance moves leftwards. This generates an 
overlap at approximately 20° (peak in Fig[S]b). These two 
regions are split again at higher inclination (Fig(8]a). Nev- 
ertheless the dynamical influence of the resonance 1/ = 31/1,2 
seems to be always predominant. 

When the semi-major axis of Saturn is greater than 
8.8AU, the dynamical situation becomes richer but more 
complicated. Indeed, the resonances of Family //generated 
by the frequency 1/1,2 are still present, but although they 
are involved in the global dynamics of the present Solar 
system (see Paper I ) , their dynamical influence decreases 
strongly. On the contrary, the secondary resonances of Fam- 
ily //associated to 1/3,7 and above all 1/2,5 begin to gener- 
ate instabilities. On one hand, the order of the 2:5 and 3:7 
MMRs being greater than the one of the 1:2 MMR, the 
corresponding resonances of Family II will probably perturb 
the Trojan swarms less than the ones depending on the 2:1 
MMR. On the other hand, one can expect to find strong 
chaos when at least two of these resonances overlap. Lets 
us notice that, because z/3,7 = 1/1,2 + 1^2,5, when a reso- 
nance of Family //associated to 1/1,2 and another one as- 
sociated to 1/2,5 overlap (intersection of red lines and green 
lines in Fig[S]d) a secondary resonance connected to 1/3,7 is 
present in the same place. But generally, the order of this 
resonance is quite high. As we have retained only the res- 
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onances v = ±1/3,7 and 2v = ±1/3,7 (the resonances asso- 
ciated to v — ±2^3,7 are also of interest, but for the sake 
of clarity they are not drawn in the figure) only four triple 
intersections are inside the studied domain (Fig.[8]d). These 
intersections are defined by the relations: 



v — 


2 


= -§^2,5 


= -f3,7 


at about 


9.04,4*7 


(16) 


v — 


2 

-Vl,2 
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= ^3,7 
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(17) 
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1 
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9.37,4(7 


(18) 


V = 


2 

-2^1,2 
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= -2^2,5 


1 

= 2^3.7 


at about 


9.5AU 


(19) 



These four different overlaps coincide with local maxima 
of the ejection rate. The most striking phenomenon occurs 
when the resonances of Family //associated to v = 1/3,7 
cross the Trojans' phase space (see formulas (|17[l and (1181) ). 
for a 2 G [9.33, 9.37] AU. Indeed, at I* = 30° (FigUa), a 
steep peak rises around 9A1AU and reaches 80%, while a 
smaller one is located on its left at 9.3AU. When the initial 
inclination falls from 30° to 2°, these two unstable struc- 
tures merge to generate the rounded peak around 9. 37AU 
(FiglSJc). The predictions given in Fig[8]d suggest that the 
left above mentioned structure is generated at 30° by the 
triple resonance defined by equation (|17|l . while the right 
one is generated by the intersection of the 2v — v\p and the 
v = —^2,5 secondary resonances. When /* decreases, the 
higher peak follows the Family //resonances associated to 
v = —^2,5 and is consequently shifted leftwards, while the 
other one is driven by the v = v\$, and goes rightwards. The 
merging of those two structures occurs at /* = 2° (Fig. [H]c) 
and generates the high rounded peak around 9.37AU (see 
eaqation HE]))- The last resonance maintains the relatively 
high level of ejection on the left side of the peak. 

In order to verify that the mechanism that has just 
been presented is the most relevant one and that it is re- 
ally Family //which drives the global dynamics of this re- 
gion, wc will present a sequence of simulations in the same 
spirit as in section 1^.2,21 Indeed, if the resonances of Family 
//play an important dynamical role, they should be visible 
in the frequency space, and particularly on the (v, g) projec- 
tion. Although we have seen in sections 13.1.11 and 13.2.21 (see 
also Paper I) that 7T( 9 , S ) T (resp. 7T(„ ]S ) o T) which is the 
composition of the frequency map with the projection on 
the (g, s) plane (resp. on the (v, s) plane) defines a isomor- 
phism from the regular regions of 23/* to 7T( s ,s) (©l* ) (resp. 
K{v,a) th e projection itu,,^ o T on the (v, g) plane is 

not one to one jGabern et a.1.1 [20051 ') . Consequently a point 
on the (y, g) plane possesses several pre-images making the 
identification of resonances very difficult. To overcome this 
difficulty, we have decided to cut the phase space by the 
plan of coordinates (a, I), limited to the rectangle 23' = 
[5.21, 5.4] x [3° , 33°] , where the other initial conditions satisfy 
( e ( ),A( Ut°)(!(»») = (e< 0) ,A< 0) +7r/3,^ 0) +7r/3,^ 0) ). 
Denoting 0' the image of 23' by the frequency map, the 
map ft(v,g\ J~ now defines a isomorphism between 23' and 
i>,s)(©')- This is clearly shown by Fig.lTOl As Fig. [5] Fig. 
UOl is composed of two blocks. The right block corresponds to 
dynamical maps of the domain 23' for eight different values of 
the initial semi-major axis of Saturn. These eight values are 



Table 4. Main resonances crossing the Trojan swarms for 02 £ 
[9.26, 9.4] AU. The first column gives the label of each simula- 
tion, these labels are reported in Fig. 1101 The initial values of 
C12 adopted in each simulation are presented in the second col- 
umn. The labels and the definitions of the resonances emphasize 
in Figs. ITOl A-H are given in the two last columns. 



igure 


a 2 (AU) 


Label 


Resonance 


A 


9.26 


1 


1/1,2 - 2u + gi = 


- 


- 


2 


i/i,2 - 2u + g = 


- 


— 


3 


1/1,2 - 2i/ - g + 2gi =0 


- 


— 


7 


Vl,3 +4i/ + 2s 2 =0 


B 


9.30421 


1, 2, 3 


defined above 






E 


1/4,9 — 3i/ + 231 + 3c/2 = 


C 


9.33368 


1, 3, e 


defined above 






8 


1/3,7 - v + 9 + 3gi =0 


D 


9.35579 


3, e, 8 


defined above 






9 


1/3,7 — v + 4#2 = 


E 


9.36316 


9 


defined above 






a 


1/2,5 + v + 2g + 32 = 


F 


9.37053 


a 


defined above 









^2,5 + v + g - 9i + 3g 2 =0 


G 


9.39263 


a, 13 


defined above 






7 


!/2,S + V + 3#2 = 






8 


^2,5 + v - g + 2gi + 2g 2 = 


H 


9.4 


P, 7, S 


defined above 


to H 


9.26 to 9.4 


Si 


s — s 2 = 






S-2 


s - S2 + 31 - 92 = 






i 


ni — 13z/ — 32 = 






5 


m — i4z/ — 3 = 






6 


ni — 14f - 59 + 4g 2 = 



given in Table [4] The color code associated to the diffusion 
rate being the same as in Fig. [5] it is not reported here. The 
left block is the corresponding view in the frequency space. 
Here again, the correspondence between the (a, /) plane and 
[y, g) plane is easy to establish. We start from L4 (bottom 
left corner in the (a, /) plane) which corresponds to the right 
vertex of 7T(i/,g)(®') where (y, g) — (8820, 362). When o 2 in- 
creases, / remaining constant, v decreases while g increases 
up to the point (v,g) = (7680, 451). On the contrary, always 
starting at L4, if / increases when a is constant, both v 
and g decrease down to the point of coordinates (8050, 270). 
The frequency map being continuous, it is easy to deduce 
where the images of the two other edges of the rectangle 
23' are located. In the frequency space (left block) vertical, 
or nearly vertical, structures are dominant. As we will see 
later, these structures are associated to resonances of Fam- 
ily II of the form v + kg + c = 0, where the absolute value 
of the integer k is lower than 2 and c is a constant real 
number. On the contrary, it is important to notice that no 
horizontal line is visible. This confirms the fact that Fam- 
ily /Vdoes not have significant influence on the dynamics 
of this region. Although there exists a priori no direct rela- 
tion between the ejection rates presented in Fig. [S] and the 
present simulations, for the reason that the first ones are de- 
duced from integrations in the (a, e) plane while the second 
ones concern the (a, /) plane, these new simulation illustrate 
very well the predictions presented in Fig. [HJd. But before 
comparing these predictions to the numerical simulations, 
let us first recall that the location of some families on either 
the (a, e) or (a, /) planes is practically independent of the 
value of the initial semi major axis of Saturn a 2 . This is the 
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Figure 10. Destabilization of the Trojan swarms by the resonances of Family II . Left block: Dynamical maps in the frequency space 
(projection on the (u, s) plane). Right block: Dynamical maps in the action space (the initial conditions of the fictitious Trojans are 
chosen in the (a, I) plane, the four other elliptic elements being fixed). The color cod associated to the diffusion rate is the same as in 

Fig.m. 



case for the secular resonances of Family III (at least in the 
studied interval [9.26, 9.4]^!7) and especially for Family I . 
Consequently these resonances are present at the same place 
along the sequence. But the intersections of the resonances of 
Family III with the (a, 7) plane are different from their inter- 
sections with the (a, e) plane. Consequently their shape will 
be different here than in section [5] While their intersections 
are similar to arches (arcs of ellipses centered in L4) in the 
(a, e) plane, in the (a, I) plane their shapes are comparable 
to the curves plotted in Fig. [4] In particular, the resonance 
s — S2 = is quite similar to the curve in Fig. [4] labeled with 
1, while the resonance s — S2 + gi — gi = is well represented 
by the curve 3. In the right block of Fig. 1101 the resonance 
s — S2 + gi — g-z = 0, denoted S2 in Fig. [10] and Table [4] , 
is always clearly visible above / = 10° while the resonance 
s — S2 — (denoted Si) stands out only in Fig. [TO] H, E 
and F. In the other plots of the sequence, Si and other reso- 
nances of Family II merge in the red region lying in the right 
part of the figures. On the other hand, the identification of 
these two resonances is not so easy in the frequency projec- 
tion on the (v, g) plane. The presence of s — S2 + gi — 32 = 
can be deduced from the existence of the curve, located at 
the bottom of the frequency domain, which corresponds to 
singularities of the frequency map. This curve is labeled with 



£?2- The resonance s — S2 = is only detectable by the in- 
stability that it generate in the (v, g) plane. But as these 
resonances do not play a significant role in the mechanism 
we want to describe, it is not necessary to give more de- 
tails. For the same reason, we will not discuss any longer 
the shape and the location of this structure. For the sake of 
completeness, three resonances of Family /are represented 
in the frequency space , labeled respectively with 4, 5 and 6 , 
for the combinations m — Viv — 32 = 0, m — 14^ — g = and 
m — I4v — 5g + 4<?2 = 0. In the (a, I) plane, the resonances 4 
and 5 are not discernible for they are located close to 5AAU 
in the same chaotic region than the resonance 52, while the 
resonance 4 is located around 02 = 5.34AU (see Fig I10I C 
right block). 

According to Fig. (8]d, if we start the migration from 
C12 = 9.26AU and increase the initial semi-major axis of 
Saturn up to 9AAU, the first noticeable resonances of 
Family //encountered are the ones defined by the relation 
1/1,2 — 2v + • • ■ — 0, where the dots represent combinations 
of secular frequencies. Indeed, in Fig. 1101 A left, three reso- 
nant lines stand out very clearly. They are emphasized by 
red lines labeled with 1, 2 and 3 (the associated combina- 
tion of frequencies is given in Table [4}. When 02 slightly 
increases, these resonances move towards higher values of v, 
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which is in perfect accordance with our predictions. The res- 
onances related to 2v = vip leave permanently the Trojan 
swarms in Fig. llOI D for a 2 = 9.35579/4(7. In the (a, I) space, 
these resonances and more generally the resonances of Fam- 
ily II appear as elliptically shaped structures centered on L4. 
In the simulations A to C, these structures move towards 
L4 before leaving the phase space. It is worth mentioning 
that these instabilities are not generated by a single reso- 
nance but by a multiplet of resonances which overlap. This 
is particularly visible in the frequency space, where parallel 
resonant lines are observable in chaotic regions (this phe- 
nomenon is detailed in Paper I ). But, for the sake of clarity, 
only one resonant line is drawn in these structures. Accord- 
ing to Fig. [8jd, when the resonances considered above leave 
the phase space (red line labeled with 1/2), the subfamily 
associated to v = 1/3,7 begins to enter the swarm. This is ex- 
actly what happens in simulation C where the black line (de- 
noted 8) corresponding to the resonance 1/3,7 — v+g+3gi = 
appears on the left hand side of the (v,g) plane (see Fig. 
UOp . As for the resonances 2v — 1/1,2 which cross the phase 
space towards L4, the subfamily associated to 1/3,7 enters the 
swarm in Fig. llOI C (label 8) and is still present in simulation 
E (label 9). Among the eight simulations that we performed, 
the simulation E, which corresponds to a 2 = 9.36316AZ7, is 
perhaps the most interesting. Indeed, just before the 1/3,7 
subfamily exits through the L4 point, the resonances de- 
noted a, j3, 7 and 5 enter the swarm. These new resonances 
corresponding to the relation v — —1/2,5 move in the oppo- 
site direction than the previous ones. This corresponds to 
the fact that, in Fig.[8]d, the green curve (labeled with — 1) 
associated to v — —1/1,2 is decreasing while the previous one 
(the black curve labeled with 1) is increasing. It turns out 
that these two structures collide giving rise to the broad 
red arch surrounding L4 (label 9 and a in Fig. IIOI E). The 
three last simulations show the evolution of the resonances 
of this last subfamily when 02 increases up to 9AAU, which 
validates one more time the predictions given in section [3T3T 
The study of this region proves definitely that the res- 
onances of Family //dominate its dynamics, and that, as it 
was mentioned in the beginning of this section, this family 
generates the main part of the peaks of instability present in 
Fig.UirC. Before completing this section, let us mention that 
two other resonances, which are rather anecdotic, appear in 
the previous simulations. Resonance 7, which is defined by 
1/1,3 + 4i/ + 2s2 = is observable in simulations A and B. 
In fact, in Paper I we have already encountered resonances 
of Family //associated to the frequency 1/1,3 at high incli- 
nation. On the other hand, it is much more surprising to 
realize that the resonance 1/3,7 — v + g + 3gi = denoted 
e stands out against the background of the frequency space 
and generates its own unstable structure visible on the (a, /) 
plane at least in Fig llOI C, 

3.4 Dependence on the initial phases 

Adiabatic invariance theory suggests that, during migration, 
the p lanets are never captured in MMRs (|Morbidelli et all 
2005). For this reason, and even if we do not take migration 
into account, we will study the evolution of the resonant 
structure along the "hyperbolic segment" defined in section 
13.1.21 Contrarily to the "elliptic segment", the former has 
the advantage to cross the 2:1 resonance at its smallest sec- 
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Figure 11. Dependence of the ejection rate on the initial phases 
for the resonances of Family II : the bottom frame shows the 
predicted locations of the resonances v 3i/i,2, v sa 2i/i,2, v ~ 
3/2i/i,2 and v sa 1/1,2 (see Fig. [8}d), for the "elliptic section" in 
solid lines and for the "hyperbolic section" in dashed lines. The 
top-frame shows the ejection rate for these two sections. 

tion, minimizing the range of Saturn's semi-major axis for 
which the two planets are trapped inside the MMR. This 
simulation is probably closer to what is expected by adding 
planetary migration. 

Denoting by 1/1,2 the quantity ni — 2n2 evaluated along 
the hyperbolic segment (Fig[2] green curve) and v^ 2 t ne one 
computed on the elliptic path (Fig[2]red curve), we deduce 
from Fig[5]that when a 2 G [8.3, 8.85] (outside of the MMR), 
the inequality < i/f^ holds. Moreover, the greater the 
value of 02, the smaller the value of — V\$,. In addition, 
the width (in a 2 ) of the 2:1 MMR along the two sections 
are very different: the "elliptic width" is about three times 
greater than the hyperbolic one (see also Fig.[T}. As a result, 
the slope of the red curve is always greater than the slope of 
the green one. These facts have direct consequences on the 
secondary resonances of Family //associated to 1/1,2. Fig llll 
illustrates the comparison between these two choices of ini- 
tial phases. In its bottom frame, the frequencies 3,2,3/2 
and 1 times 1/1,2 are represented by solid lines for 1/^2 and 
by dashed lines for i/^ 2 - Just like Fi^5}d, this plot enables 
us to predict the locations of the resonances of Family II . 
The associated numerical simulations are reported in the 
top frame of Fig llll FI. the ejection rates along the hyper- 
bolic section are represented by a dashed curve with black 
circles while the solid line with white circles corresponds to 
the elliptic section. Again, the numerical simulations and 
the predictions are in very good agreement. As regards the 
location of the resonances of Family II , a given resonance is 
crossed by the Trojan swarm at a smaller value of a 2 in the 



' In this section, the simulations arc limited to the Trojans having 
an initial inclination equal to 2°, since one example is sufficient 
for the purpose of studying the dependence of the results on the 
planets' initial phase. 
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hyperbolic case than in the elliptic one, and the splitting 
between these two locations decreases when the planetary 
system moves away from the orbital resonance. The width 
of the resonances of Family II is also affected by the change 
in the initial phases in the planets. Because the slope of 
the curves plotted in the bottom frame of Fig. [11] is larger 
for the elliptic section than for the hyperbolic section, the 
resonances generated in the hyperbolic case are wider (i.e. 
they occur within larger intervals Aa) than the ones asso- 
ciated to the elliptic case. But, as shown in Fig. 111! even 
if the ejection peaks are wider and shifted leftwards along 
the hyperbolic section, the global features are the same in 
both cases. Consequently, the modification of the planetary 
phases does not introduce appreciable effects on the global 
dynamics of the Trojan swarms. 



4 DISCUSSION 

The global dynamics of a Trojan swarm is shaped by its res- 
onant structure. This structure depends crucially on the fun- 
damental frequencies of the planetary system in which it is 
embedded. By generalizing the results of Paper I concerning 
Jupiter's Trojans, we show that in a general planetary con- 
figuration this resonant structure can be decomposed in four 
different families. Among these families, two are particularly 
sensitive to the geometry of the planetary system: Family 
II and Family IV axe dependent on the closeness of the plan- 
etary system to MMRs. It turns out that variations of the 
semi-major axis of the planets deeply modify the dynamical 
influence of these two families. 

Based on this decomposition and on the dependence 
of the families on fundamental frequencies of the planetary 
system, we present a general method making possible the 
study of the global stability of Trojan swarms during plan- 
etary migration possible. We show that the knowledge of 
the evolution of the fundamental frequencies of the system 
during its migration suffices to predict the main features of 
its global dynamical evolution (transition between total sta- 
bility and strong instability), and to find which planetary 
configurations which lead to stable or unstable swarms. 

In our application to the case of the Jupiter Trojans 
disturbed by Saturn, we have varied one parameter: the 
semi-major axis of Saturifl While in a realistic migration, 
many other parameters would vary, this would modified only 
slightly the results obtained in this paper. For example, the 
variations of the planetary phases have almost no effect 
on our results (section I3.4[) . In the same way, the "apsi- 
d al corotation" of the plan etary perihelia which, according 
to lMarzari fc Scholll (I2007T ). strengthens the influence of the 
secular resonances, probably does not significantly modify 
the resonant structure of the Trojans. In section 13.2.21 we 
have described the mechanism leading Family IV to generate 
strong instabilities when the planetary system approaches a 
MMR. 

Section 13.31 analyzes the dominant role of Family II in 



Equivalent results are obtained if one varies, instead, both 
semi-major axes, i.e. the results depend only on the ratio of the 
two axes. 



the production of unstable regions during planetary migra- 
tion. In particular, we demonstrate a mechanism of forma- 
tion of huge chaotic regions by the merging of several reso- 
nances crossing the Trojan swarm in opposite directions. 

The resonant structure, and particularly its splitting in 
four families as defined in section 12.31 is generic. Moreover, 
the definitions of the families seem to hold for almost all 
planetary configurations. But the respective dynamical 
influences of the four families depend on the mass of the 
planet harboring the Trojan swarms. According to formula 
{TJ, if this mass decreases, the decrease in the proper 
frequency g is proportional to the planetary mass while the 
one in v is weaker (proportional to the square root of the 
planetary mass). The elements of Family IV being approxi- 
mated by the relation: g ~ v v<q , the smaller the planetary 
mass, the closer this resonance is to the p : q MMR. As a 
result, the smaller the mass of the planet, the more local 
the dynamical influence of Family IV . In this case, we can 
expect that the role of the resonances of Family IV becomes 
negligible compared with the role of Family II . Conversely, 
the secular resonances (Family III) which do not involve 
the frequency s, seem to be negligible in the Jovian Trojan 
swarm. But their contributions increase for smaller plane- 
tary masses. This has been a lready men t ioned in section [231 
This has been observed by iBodossianl (|2008h (available at 
http: / / www.imcce.fr /page.php?nav=en /publications / theses / 
index. php) in the case of the Trojans of Saturn. In both 
cases, Family II and Family IV dominate the dynamics of 
the swarms during migration. 

Finally, the method of prediction developed in this pa- 
per being very general, it can be used to study the dynamics 
of a large class of problems. Beside its straightforward ap- 
plication to the Trojans of a planet (results regarding the 
dynamics of the hypothetical Trojans of Saturn will be pre- 
sented in a forthcoming paper) , it can also be used to study 
the dynamics of satellites in a planetary system in migration, 
or more generally, to understand the behavior of a dynamical 
system undergoing quasiperiodic perturbations with slowly 
varying frequencies. 
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